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INTRODUCTION 


En paléontologie, les exigences actuelles de la détermination spécifique 
sont celles d’une discrimination toujours plus précise et plus fine, exercée, non 
plus sur des individus isolés, mais sur des collections abondantes. Or, à la base 
de ce travail les éléments originaux des genres ou espèces choisis par de nom- 
breux chercheurs sont souvent d’une incroyable déficience. Ces types ou holo- 
types ont, pour la plupart, été sélectionnés très anciennement, souvent d’après 
le plus bel échantillon et dans l'ignorance complète des formes moyennes corres- 
pondant à chaque phase de la croissance, ainsi que des formes déviées par l’action 
du biotope. Aussi de nombreux spécimens types sont-ils devenus insuffisants 
pour les études comparatives et descriptives. De plus beaucoup d’entre eux, 
éminemment précieux pourtant, ont été dispersés, soit au cours de leur pérégri- 


nation dans les institutions, soit par la négligence des collectionneurs. 


À côté de ces indigences, dues au facteur humain, il en est d’autres, inhé- 


rentes à l’objet lui-même-et dues à son polymorphisme et à sa fossilisation. 


Les Pélécypodes des faunes limniques ou dulcicoles des terrains houillers 
offrent un exemple magistral de la faible valeur de leurs types utilisés jusqu'à 
présent. Ils ont été prématurément choisis au sein des collections numérique- 
ment très insuffisantes et non véritablement représentatives des espèces. Cette 
situation a créé de multiples confusions dans les travaux paléontologiques par 
l'établissement de diagnoses hâtives et de déterminations approchées, qui allon- 
gent et remanient sans cesse la liste des synonymies. Il faut ajouter qu'une icono- 
oraphie défectueuse complique encore la tâche du systématicien. 
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Ces coquilles de Pélécypodes sont, de près, toutes sembläbles et, de très 
près, toutes différentes par la variabilité de leurs formes résultant des actions 
du milieu. Leur étude comparative par le simple examen visuel, ayant pour but 
leur classement en groupes, soit leur description morphologique, présente des 
insuffisances et est, au surplus, fallacieuse. Aussi les nécessités de se faire bien 
comprendre ont-elles amené déjà quelques chercheurs à traduire les résultats de 
leurs observations par des mesures. 


Allant plus loin, on peut penser que l’abstraction mathématique pourrait 
fournir aux formes envisagées un cadre, des limites — puisqu'il n'apparaît pas 
de limites biologiques reconnaissables — qui soient à la fois larges et précises 
et suffisamment claires pour contenir l'espèce et ses variétés. Dans l’utilisation 
des nombres à la définition de l’espèce, rares seraient sans doute les cas où l’on 
pourrait envisager la « géométrisation » de sa forme et il est certain que celle-ci 
ne devrait être appliquée que là où les bases morphologiques de la détermination 


conchyliologique sont peu nombreuses, simples et mesurables avec précision. 


Cette voie apparaît, dans l’avenir, riche de possibilités pour le classement 
systématique de ces fossiles, quand on envisage les travaux, récents ou en cours, 
sur la définition des coquilles par l'établissement de leur équation surfacielle 
basée sur la courbure du bord libre [L. Lison (°)|. 


Au préalable et présentement il est nécessaire d’accroître d’abord considéra- 
blement les collections de fossiles des espèces limniques ou dulcicoles de 
_l'Anthracolithique; ensuite, d'opérer un classement provisoire et très prudent 
de ces formes sur les bases des caractères morphologiques simples et des carac- 
tères régionaux et stratigraphiques. C’est dans ce sens qu'a été conçue cette 
« Étude biométrique des Anthraconauta du houiller de la Belgique », dont la 
première partie est consacrée à l'examen comparatif d’Anthraconauta bel- 
gica Hip et d’'Anthraconauta minima (Hinn, non Lupwic), sur les bases morpho- 
logiques usuelles, mais à l’aide de méthodes mathématiques et statistiques bien 
définies. 


() LisoN, L., 1942. 
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Étude biométrique des Anthraconauta du Houiller 
de la Belgique 


(PREMIÈRE PARTIE) 


CHAPITRE I. 


$S 1. Place des mollusques limniques ou dulcicoles du Terrain houiller 
dans la classification zoologique systématique. 


Il n’est pas sans intérêt, dans une étude sur les Anthraconauta, de passer en 
revue l'opinion des malacologistes et des paléontologistes sur la place occupée par 
les mollusques limniques ou dulcicoles de l’'Anthracolithique dans l’une ou l’autre 
classification. Très peu nombreux sont ceux qui, depuis W. Hip, ont fait œuvre 
de systématicien en étudiant de près ce problème, à propos duquel on trouve plus 
de redites que de recherches originales. 

Trois courants d'opinion ont cours. Le premier classe les genres Anthra- 
comya, Carbonicola et Naiadites dans la famille des Unionidae. Le second, qui a 
rallié de nombreux avis, rapporte ces genres, ou les deux premiers de ceux-c1, 
aux Cardiniidae. Le troisième range ces genres dans une famille propre, tout en 
reconnaissant la diversité des affinités qui y règne. 

M. Neumayr (°) (*), dont la classification est basée sur les caractères de la 
coquille, principalement ceux de la charnière, a rangé Anthracomya SALTER dans 
les Unionidae et rapproche ces dernières des Trigonidae, alors que H. Dou- 
VILLE (*), venant après lui, faisait des Anthracosia des Unionidae vraies, issues de 
Taxodontes. Pour W. Hip (°), Anthracomya Sarrer et Carbonicola M'Coy sont à 


(2) (5) Neumayr, M., 1893 et 1891, pp. 412 et 719. 
(*) Douvicee, H., 1912, p. 441. 
(5) Hip, W., 1894, pp. 12-38 et suiv. 
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placer dans la famille Unionidae, tout en notant que les Anthracomya possèdent 
un ensemble de caractères qui les rapprochent plus des Mytilidae que les Carbo- 
nicola, lesquelles seraient intermédiaires entre les premières et les Unios 
modernes. Quant aux Naiadites, rapportées tour à tour par les auteurs aux Modiola, 
aux Mytilus, aux Dreissensia (par Lupwie) et aux Myalina (par BarRotrs), Hip les 
rapporte aux Mytilidae. 

DE KoniNG et ne RYcKkHoLT ont considéré les « Unios du Houiller » comme 
affines du genre Cardinia d'Acassiz. K, ZiTTEL, G. STEINMAN et L. DoperRLEIN (°) 
ont rangé Anthracosia Kixc dans la famille des Cardiniidae, ordre des Hétéro- 
dontes. Dans la classification de W. H. Dai (°), on{rouve encore Carbonicola 
M'Coy, Anthracosia Kinc et Naiadites DAwson rassemblés dans la famille des Car- 
diniidae de ZITTEL, superfamille des Naïadaceae, ordre des Schizodontes. C'est à 
cette classification que se sont ralliés la plupart des paléontologues anglo-saxons, 
entre autres J. HEnpErson (°), et, en Belgique, M. F. DEMaAnerT (°). P. PELsE- 
NEER (°°), dans sa classification zoologique basée sur les branchies, a maintenu le 
genre Anthracosia Kinc dans la famille des Cardiniidae de Z1rTTEL, mais a rangé 
celle-ci dans le sous-ordre des Submytilacea, ordre des Eulamellibranches. Il a 
classé le genre Anthracoptera SALTER dans la famille des Myalinidae, sous-ordre 
des Pectinaceae, qu'il distinguait des Mytilaceae, ces deux derniers sous-ordres 
faisant partie des Filibranchiés. F. BernarD (7) a constaté, sans prendre position, 
que les Anthracosia Ki et les Carbonicola M’Coy sont rapportées tantôt aux 
Cardiniidae, tantôt aux Unionidae, mais que, dans l’un ou l’autre cas, elles 


seraient des formes ancestrales communes aux deux groupes, comme le signa- 
laient PonLre et ZITTEL. 


Parmi les chercheurs partisans de réunir les Lamellibranches limniques ou 
dulcicoles du Houiller en une famille propre, il faut citer W. Amazizxy (©), qui 
groupa dans la famille des Anthracosiidae (ou Anthracosinae) les genres Carbo- 
nicola M'Coy, Anthracosia Kinc, Palaeomutela Amarizxy, Oligodon Amarrzxy et 
Naiïadites Dawson. Pourtant, il reconnaissait que les deux premiers genres sont 
des formes hétérodontes et que le genre Naiadites est une forme dépourvue de 
charnière. Pour M. P. Pruvosr (*), Carbonicola M'Coyx et Anthracomya SaLrer 
doivent être rangés dans la famille des Anthracosiidae d’Amarizxy, et Naiadites 


(5) STEINMAN, G. et DODERLEIN, L., 1890, p. 265. 

(7) DaLz, W. H. in Z1TTEL, K. and EASTMAN, R., 1927, pp. 451-452. 
(5) HENDERSON, J., 1935, p. 291. 

(°) DEMANET, F., 1938, pp. 139-145. 

(1°) PELSENEER, P., 1906, pp. 262-267. 

) BERNARD, F., 1895, p. 581. 

12) AMALIZKY, W., 1892, pp. 136-142. 

13) PRUVOST, P., 1930, pp. 227-247. 
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Dawson et Anthraconauta Pruvosr, dans la famille des Mytilidae de LamarcKk. 
-L'homogénéité de la famille des Anthracosiidae a été reconnue encore par 
H: Pau (”). Une opinion récente de M. A. E. TruEman place les Carbonicola du 
groupe de C. aquilina Sowergy dans la famille des Anthracosiidae (‘). 


$S 2. Phyllogénèse des Anthraconauta. 


Les opinions émises ces dernières années sur l’origine phyllétique des 
Anthraconauta traduisent une grande diversité de conception. M. D. Feporov (°°) 
rattache les grandes Anthraconauta du bassin de Kusnetzk aux Anthracomya à 
faible obliquité. De son côté, M. H. Wenri (7) a exprimé l’idée d’une descen- 
dance phyllétique issue des Naiadites (*), opinion qui a subi la critique de 
MM. E. Dix et À. E. TruEman (°°). Ces deux thèses ont été discutées récemment 
par M. J. Werr (°°). Ce paléontologue a signalé la possibilité de trouver une 
origine phyllétique des Anthraconauta dans le groupe des Posidonia de la famille 
des Pteriidae (Aviculidae). Il estime cependant possible, à l’origine, une sépara- 
tion diphyllétique des Anthraconauta et des Naïadites (provenant les unes et les 
autres d’un stock commun au stade « Modiola ») ÉD PRnTE différenciées et 
stabilisées en une morphologie distincte (*”. 

Cette suggestion avait été faite auparavant par M. N. D. NEWELL (*), qui a 
étudié d'une façon approfondie le problème des relations entre Myalinidae et 
Mytilidae du Paléozoïque de l'Amérique du Nord. Après avoir constaté que la 
tendance à l’obliquité des valves et le déplacement antérieur de l’umbo n'est pas 
une tendance particulière aux Myalinidae, mais se retrouve également chez les 
Mytilidae, les Ambonichiidae et les Pteriidae (*), il conclut à un cas d’homéomor- 
phisme, non nécessairement phylogénétique, entre les formes de ces groupes. De 
plus, il reconnaît que les Mytilidae primitives du Paléozoïque présentent une 
expression « modioloïde » très nette, ce qui l’incite à classer la majorité de leurs 
formes ayant un aspect « mytiloïde » ou « modiloïde » parmi les Mytilacea, sous 
quelques réserves cependant. | 

Il lui apparaît que; par la forme et la structure de la charnière, le genre 
(4) PAUL, H.. 1941, pp. 276-277. | 
(5) TRUEMAN, A. E., 1941, communication écrite; voir aussi : TRUEMAN, A. E., 1946. 
(#5) FEpoTov, D., 1937. 

(7) WEBRLI, H., 1933, pp. 133-134. 

-(#) WEIR, J., 1944, pp. 316 et suiv. 

() Dix, E. et TRUEMAN, À. E., 1931, pp. 192-195. 
(2°) WEIR, J., 1944, p. 36. 

(1) Wei, J., 1944, p. 321. 

(2?) NEWELL, N. D., 1940. 

(2) NEWELL, N. D., 1940, p. 8. 
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Modiolopsis HALL, déjà présent à l'Ordovicien, pourrait êtré vraisemblablement 
l'ancêtre de tous les Mytilidae. Les espèces de ce genre sont proches des Volsella 
(Modiolus) actuelles (?*). 

Cette opinion, à laquelle nous nous rallions, reprend des idées émises ancien- 
nement par F. BERNARD (”), qui voyait dans les Modiola et les Modiolopsis 
d'ULricu des époques précarbonifériennes, les stades divers du développement 
des Mytilidae. 


8 3. Caractères généraux des Anthraconauta. 


C'est à M. P. Pruvosr que revient la création du genre Anthraconauta, qu'il 
rattache à la famille des Mytilidae de Lamarcx en séparant les espèces de ce genre 
du groupe générique désigné sous le nom d’Anthracomya, famille des Anthra- 
cosiidae d’AmaLizxy. Cet auteur a exposé longuement les raisons qui l'ont amené 
à établir cette coupure, dans un travail sur la faune continentale du Westphalien 
de la Belgique (*°), dont voici extraite la diagnose du genre : « Valves de contour 
subovale, allongées obliquement (angle d’environ 80°) par rapport au bord 
cardinal, toujours plus court que la plus grande longueur. Large échancrure 
byssale peu profonde sur le bord antéro-ventral, séparant la valve en une partie 
antérieure arrondie en oreille et une postérieure de beaucoup plus importante. 
Une carène mousse rectiligne, partant du crochet, forme la ligne de plus grande 
longueur. Caractères internes inconnus. Test de structure prismatique, recouvert 
d'un périostracum épais, souvent ridé, portant des stries d’accroissement concen- 
triques irrégulièrement espacées. Formes vivant associées, en colonies, flottant 
fixées par leur byssus à un support » (?7. 

À l’époque de la publication de l’ouvrage, l’auteur avait reconnu l'existence 
des cinq espèces suivantes : A. Phillipsi (WiLriamsoN), génotype, À. minima 
(Lupwic), À. laevis (Dawson), A. belgica (Hinp) et A. stephaniensis Pruvosr. 
Depuis, la liste s’est augmentée de : A. tenuis Davies et TRruEMaAN (*), 4. Wrighty 
Dix et TRUEMAN (*), A. anthracomyoides Fenorov (°°), À. krochalovskiensis 
Fenorov (*”), À. tenuistriata Fenorov (*?), À. iljinskiensis Feporov (**), A. pseudo- 
phillipsi Feporov (9, À. kemeroviensis Feporov (**), À. gigantea Racozin (°°), 


(2) NEWELL, N. D., 1940, pp. 35-36. 

(*) BERNARD, F., 1898, pp. 196-197. 

(*) PruvosT, P., 1930, pp. 223-225. 

CD RPEUVOST 2 MOS0 ED 27e 

(25) DAVIES, J. H. and TRUEMAN, A. E., 1927. p. 235. 
(2°) Dix, E. and TRUEMAN, A. E., 1931, p. 184. 

(80) (81) (#2) FEDOTOV, D. M., 1938, pp. 219-250. 

(35) (54) (55) FEDOTOV, D. M. 1937, pp. 103-104. 

(55) RAGOZIN, U., 1933, pp. 56-63. 
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A: sibirica Racozi (°”), A. minuta Racozin (°°), A. longa Racozn (°°), A. flagi- 


« 40 . à 
_tiosa Racozin (”), A. pungabica Reep (”), toutes espèces peu connues encore. 


Récemment, le terrain houiller du Lancashire a permis à M. Dewar de distin- 
guer les espèces et variétés suivantes : À. minima (Hivp, non Lupwic), var. 
obliqua, À. trapeziforma Dewar, A. trapeziforma Dewar var. tenuides Dewar, 
A. subovata Dewar, A. subovata var. candela DEWAR, et de redécrire À. minima 
(Hip, non Lupwic) et À. scotica Erneripce (*). 


Il faut noter que quelques auteurs, H. Weurir et P. Kuxux, continuent de 
ranger les Anthraconauta distinguées par M. P. Pruvosr parmi les Anthracomya. 


8 4. Spécificité d’'Anfhraconauta belgica Hinp. 


Parmi les espèces distinguées par M. P. Pruvosr, l’Anthraconauta belgica 
Hip est une des moins connues, parce que de rencontre très rare. Tous les 
gisements qui en ont livré des exemplaires sont situés dans le bassin houiller de 
Liége, en Belgique. L'examen des quelques échantillons que possèdent les collec- 
tions de nos institutions a éveillé des doutes sérieux sur la spécificité de l’espèce, 
ce qui nous a conduits à entreprendre une revision critique du sujet. 

Voici la diagnose originelle établie par W. Hip : « Coquille petite, ovale- 
ment subcirculaire; l'extrémité antérieure est courte, profonde, avec un bord 
régulièrement courbé. Le bord inférieur est presque droit et descend rapidement. 
L'angle postérieur, inférieur, est très dilaté et brusquement arrondi. Le bord 
postérieur, court, presque droit, est tronqué obliquement à partir du haut vers 
le bas. La charnière est droite, courte et postérieurement élevée. L’angle posté- 
rieur, inférieur, est très dilaté et brusquement arrondi. Le bord postérieur, 
court, presque droit, est tronqué obliquement à partir du haut vers le bas. La 
charnière est droite, courte et postérieurement élevée. L’angle postérieur est bien 
marqué et obtus. Les crochets, petits et vaguement marqués (inconspicuous), 
sont placés dans le tiers antérieur de la valve. La coquille est très dilatée et 
aplatie postérieurement avec une dépression marquée le long de la pente dorsale. 
La valve est ornée de stries d’accroissement concentriques, fines et serrées. 
Périostracum ridé » (“*). Il ressort de cette description que les caractères princi- 
paux de l'espèce seraient : la forme orbiculaire plus arrondie que celle des autres 
Anthraconauta, la région postérieure moins développée et le bord dorsal moins 
oblique. D’après les.observations de M. P. Pruvosr, « la position subcentrale de 


(57) RAGOZIN, U., 1933, pp. 56-63. 

(28) (58) (4) RaGoziN, U., 1935, pp. 57-74. 
(4) REED, F. R. G., 1939. 

(42) Dewar, W., 1939, pp. 54-61. 
CSRENINDEANVe 1011 p;t17. 
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son crochet et par suite le développement de l'aile antérieure en font un type 
distinct des autres » (). Ces caractères spécifiques ont également été reconnus 
par M. F. Demaner () (9. 


8 5. Nécessité d’une étude biométrique. 


Un examen attentif et soigné des spécimens types décrits et figurés par 
W. Hp (7) et reproduits ailleurs par M. F. DEMANET (*°) (”) révèle ce qui suit : 

Le spécimen des figures Let la, reproduites par W. Hinp, est le moulage 
d'une empreinte incomplète dans la région de l’angle postéro-inférieur, où il 
manque 1 à 2 mm. pour présenter un contour valvaire normal, parallèle aux 
stries concentriques d’accroissement. De plus, dans la région cardinale, le test 
a été plié et retroussé dans la plaquette schisteuse, ce qui a eu pour conséquence 
une déformation de la ligne de la charnière ne permettant pas une bonne appré- 
ciation de la position de l’umbo. La valve est fortement aplatie, fait qui contribue 
à lui donner un aspect plus arrondi, 


Le spécimen des figures 2 et 2a, du même auteur, a subi aussi un aplatisse- 
ment déformateur du contour valvaire. Au demeurant, le meilleur de ceux repro- 
duits par ces chercheurs est l’exemplaire de la figure 3. Il se présente, en relief, 
non déformé, avec un contour valvaire net et bien délimité sur les deux tiers de 
son pourtour; seule la région du crochet est imprécise. C’est sans conteste une 
empreinte d'Anthraconauta. 


D'autres spécimens, catalogués sous le nom d’A. belgica, les uns récoltés 
par M. F. DEMANET aux gisements d'Angleur et d'Andenne; les autres, des 
diverses provenances énumérées au $ 6, ont été examinés. Une étude soignée de 
tous ces exemplaires nous a incités à mettre en doute la spécificité de À. belgica, 
parce que le caractère distinctif principal, l’orbicularité des valves, apparaît le 
plus souvent n'être qu’accidentel, comme un effet de leur écrasement. Les exem- 
plaires rangés sous ce nom ne nous ont pas paru, à l'examen visuel, différents 
de ceux définis comme A. minima. Une étude biométrique, entreprise sur un 
matériel beaucoup plus abondant que celui récolté jusqu'alors, s’imposait pour 
confirmer ou infirmer cette observation. 


(IPPRUVOST, P:21930/ DP227. 
(#5) DEMANET, F., 1938, p. 144. 
(#5) DEMANET, F., 1941, p. 261. 
(47) Hip, W., 1911, pl. I, fig. 1-3. 
() 
(®) 


#8) DEMANET, F., 1938, pl. 119. fig. 11 à 13. 
4) DEMANET, F., 1941, pl. 15, fig. 38. 39. 
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S 6. Position stratigraphique des gisements et provenances 
du matériel étudié. 


Les gisements connus appartiennent tous aux terrains des âges namu- 
riens (Nm) et westphalien inférieur (Wnza). Voici les localités d'origine des 
exemplaires examinés et qui représentent pratiquement l’ensemble des spéci- 
mens existant dans les collections de Belgique : 


ÉTAGE NAMURIEN. — ASSISE D'ANDENNE (Nm2). 

PI. Chenée : 

1. Angleur, en face de la fabrique de produits chimiques. [Coll. Mus. roy. 
Hist. nat. HU. Bruxelles, coll. X. Srainier, 1891-1901, I.G. 11.312, dét. 
X. SraniEr : Modiola carlottae ROEMER; rév. F. DEMANET, 1939 : A. belgica 
(Hino).| 

2. Kinkempois (Angleur), gisement non précisé. (Coll. Université de Liége, 
coll géologique, n° 92, dét. X... : Mytilus carlottae RoEMER.) 

5. Angleur, gisement non précisé. [Coll. Mus. roy. Hist. nat. Belg., 
Bruxelles, anc. coll. de Konnex et Ruror. I.G. 8573. Coll. des spécimens figurés, 
dét. W: Hinp, 1902 : Anthracomya belgica sp. nov.; rév. P. Pruvosr, 1921 : 
Modiola belgica Hinp; rév. F. DemanerT, 1930 : A. belgica (Hip). Plusieurs 
spécimens types, reproduits in : Hip, W., 1911, t. VI, fig. 1-3, et in : DEma- 
NET, F., etc., 1938, pl. 119, fig. 11-13.] 

4. Angleur, à flanc de colline, à 500 m. au S.-0. de la Station. Points 
n® 21.506 et 21.508. (Coll. Mus. roy. Hist., nat. Belg., Bruxelles, I.G. 12.527, 
coll. des spécimens figurés, dét. K. DEMANET : À. belgica Hip, spécimen type 
reproduit in : Demaner, F., 1941, pl. XV, fig. 38 et 39.) 

5. Angleur, à flanc de colline, à 180 m. au S.-E. du tunnel-abri de la 
rue Jean Volders. (Coll. Ass. Étude Pal. Strat. Houil., I.G. 39, 68.) 


PI. Andenne : 

Andenne, bois de Stud, carrière. Point n° 19. (Coll. Mus. roy. Hist. nat. 
Belg., 1.G. 11.516, dét. F. Demaner, 1939 : A. belgica Hinn. Coll. Ass. Étude 
Pal. Strat. MOULE LG. n° 54, 66.) 


P1.. Herve : 

Val-Dieu, ferme de la Mouderie. Gisements n° VIII et XVII. [Coll. Ass. 
Étude Pal. Strat. Houil., I.G. n° 70, dét. Cn. Ancron et W. Van LECKWIICK : 
A. belgica (Hin»).| 


Charbonnage de Forte-Taille. Sondage n° 26 (4) à Jamioulx. A 244, 246, 247 
et 248,50 m. de profondeur, niveaux n°” 11 et 12 de la coupe. [Coll. Mus. roy. 
Hist. nat. Belg., Bruxelles, coll. X. Srainier, [.G. n° 11.312, dét. X... : Posido- 
niella sp., rév. F. Demaner, 1939 : A. belgica (Hinn).] 
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ÉTAGE WESTPHALIEN. — ASsisE DE CHÂTELET. (Wnia, Wn1b). 


Charbonnage de Courcelles-Nord, puits n° 8, sondage intérieur, de 861 à 
363 m. de profondeur, niveau n° 13, n° 157 de la coupe. [Coll. Mus. roy. Hist. 
nat. Belg., I.G. 14.632, dét. P. Pruvosrt, 1930 : A. belgica (Hin»).] 


Charbonnage d'Amercœur, siège : Naye-à-Bois, étage de 790 m., à 96,80 m. 
au Nord du puits, à 50 m. sur veine Sainte-Barbe de Floriffoux. [Coll. Mus. roy. 
Hist. nat. Belg., [.G. 14.632, dét. P. Pruvosr, 1930 : À. belgica (Hinp).] 


CHAPITRE II. 


INTRODUCTION. 


Les grandes lignes d’une étude biométrique de populations d'individus de 


provenances différentes et soupçonnés appartenir à une même espèce sont les 
suivantes : 


s 


A. — Choix des caractères morphologiques à mesurer. Grandeur de la pré- 
cision des mesures. 


B. — Choix des méthodes de la statistique mathématique à appliquer à 
l'étude biométrique. 


C. — Mensuration des individus de chaque population. 


Pour chaque population : 
D. — Étude de la corrélation des caractères. 


E. — Étude des indices caractéristiques (rapports entre les mesures des 
caractères des individus) et de leur distribution de fréquences. 


F. — Conclusions biologiques relatives à la « variabilité individuelle » et 
aux « variations spécifiques ». 


G. — Comparaison des indices types des différentes populations pour la 
recherche d’une conspécificité éventuelle de leurs éléments. 


Dans le présent chapitre, nous exposerons les trois premiers points du 
schéma précité en développant particulièrement l'exposé de la méthode 
employée. Il ne convient pas ici de s'étendre sur toutes les méthodes mathéma- 
tiques qui sont susceptibles d'être utilisées en biométrie. Signalons toutefois la 
grande indigence dont font preuve nombre de biologistes qui n’utilisent pas 
pour leurs recherches toutes les ressources de l’analyse mathématique et de la 
statistique, ou qui mettent en pratique, sans discernement, des recettes toutes 
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faites, telle la célèbre loi normale de Gauss qui n’est au demeurant qu'un 


cas particulier de la distribution des fréquences. Aussi développerons-nous les 


notions mathématiques utilisées dans cette étude pour qu'aucun doute ne subsiste 
ni sur l'opportunité d'adopter certaines méthodes plutôt que des formules toutes 
faites, ni sur la nécessité de pousser dans toute sa rigueur et dans le détail l’ana- 
lyse statistique et mathématique des séries. 


A. — CHOIX DES CARACTERES MORPHOLOGIQUES À MESURER. 
: GRANDEUR DE LA PRECISION DES MESURES. 


CARACTÈRE GÉNÉRAL DES MENSURATIONS EFFECTUÉES. — Les problèmes de 
caractérisation et d'ajustement de séries statistiques que s'efforce de résoudre la 
statistique mathématique, et que nous utiliserons, portent sur des variables 
dénommées : variables statistiques. Ces variables ne sont pas autré chose que 
l'expression quantitative des caractères observés et mesurés, Ces caractères sont 
susceptibles de varier dans certaines limites de deux façons : ou bien la variable 
statistique ne peut prendre qu'un nombre fini de valeurs formant une suite 
discrète ou discontinue : æ,, Æ, …, æ,; ou bien la variable est une variable 
continue, pouvant prendre toutes les valeurs comprises entre les limites extrèmes 
de la variation. 

Dans beaucoup de cas un grand nombre de caractères biométriques, et 
notamment les tailles d'individus, les mesures de longueurs et d’angles caracté- 
ristiques des coquilles de lamellibranches, sont des variables statistiques con- 
tinues. 


GRANDEURS MESURÉES. — La figure 1 représente les grandeurs mesurées sur 
les valves du type de l’Anthraconauta et leur notation par MM. P. Pruvosr, 
J. H. Davss et A. E. TRuEMAN et par nous-mêmes. En conclusion de cette étude, 
nous montrerons que certaines de ces grandeurs ne sont point caractéristiques 
et sont dépourvues de signification précise. 

Pour P. Pruvosr (*”), dans la figure 1 (1), les grandeurs : 

c—longueur de la ligne cardinale, 

b=longueur de la carène, 

a-b=la plus grande largeur. 


Pour J. H. Davies et A. E. Trueman (°), dans la figure 1 (2), les grandeurs : 


L=]longueur de la valve, 
H= hauteur, 
A=—longueur de la région antérieure. 


(5) Pruvosr, P., 4920, fig. 1 à 3, p. 7 et fig. 8-12, p. 30. 
(51) DAVIES, J. H. et TRUEMAN, A. E., 1927, p. 242 
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Au cours de leur emploi dans une analyse biométrique, les mesures H, A 
et T (épaisseur) sont souvent exprimées comme pourcentage par rapport à la 


longueur, et la caractérisation géométrique de la coquille s'établit ainsi : —L; 
Else À 
Li), Lif)lc 


Pour notre étude, nous avons représenté, à la figure 1 (3), les grandeurs 
ayant fait l’objet de mesures; elles ne nécessitent pas de commentaires; nous 
avons appelé g la somme e+c. 


Pour M. Gziserr (*), les grandeurs a et b, de la figure 1 (3), sont les dia- 


a : se ONE 
mètres umbono-ventral-et antéro-postérieur, et le rapport ; 1 indice des diamètres. 


L 


0 
{ 
' 
! 
! 
! 
ll 
0 
vf 


La mesure des grandeurs précitées fut effectuée, sur projection optique 
agrandie de l’image, au moyen d’un pied à coulisse de précision et d’un gonio- 
mètre d'application. La précision des mesures linéaires fut de +0,1 mm. et la 
précision des mesures angulaires de 1 à 5°, suivant l’état du spécimen. 


B. — CHOIX DES METHODES DE LA STATISTIQUE MATHEMATIQUE 
A APPLIQUER A L’ETUDE BIOMETRIQUE (*). 


SÉRIES STATISTIQUES A UNE VARIABLE. 
$ 1. Distribution de fréquences. Représentation géométrique. 


Soit æ une variable statistique discontinue susceptible de prendre n valeurs 
distinctes x,, x, .…, æ,. Le plus souvent, en pratique, ces valeurs forment une 


suite en progression arithmétique, les différences entre valeurs successives étant 
toutes égales entre elles. 


(5?) GLIBERT, M., 1943, p. 5. 
(*) Ouvrages généraux à consulter : 


VULE, G. U. et KENDALL, M., 1945; RISSER, R. et TRAYNARD, GC. E., 1933; RISSER, R., 


1932; LINDER, A., 1945; VAN DEUREN, P., 1934-1935; FELDMAN, W. M., 1935; DARMOIS, G., 
1941. 
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Un certain nombre N>n d'expériences montrent que ces n valeurs sont : 
observées respectivement @(x:), o(x), .…, o(x,) fois. On a évidemment 


p (2) + p (de) + + p (a) — » o(æi) = N. 


Les quantités o(x;) e 


t Le sont les fréquences absolues et les fréquences 


relatives ou probabilités statistiques des x. Le nombre N sera toujours supposé 
assez grand pour qu'on puisse le qualifier de « grand nombre » au sens général 


4 
= ph) 


O LXy Lo TL; NEO EE 


qu on accorde à cette expression dans le vocabulaire des probabilités. Nous dirons 
qu'un grand nombre est un nombre tel que les conclusions de nombreuses pro- 
positions du calcul des probabilités, à caractère asymptotique, relatives à des 
nombres croissant au delà de toute limite, soient applicables de façon légitime, 
avec une erreur négligeable, mais bien connue. 


On obtient ainsi une distribution des fréquences observées suivant les valeurs 
de la variable statistique. Cette distribution est susceptible d’être représentée 
géométriquement, sous forme d’un graphique à n points isolés rapportés à un 
système d’axes x et y—=o(x). Ces points P,, P,, .. P;,, P, ont pour coordonnées 
æ, et o(x,). Relier ces points par une courbe continue ou par une ligne brisée n’a 
pas de signification analytique, la fonction (x) n'étant définie qu'en #,, æ,, .…. æ, 
Mig. 2). 1: 

Considérons à présent une variable statistique continue x, susceptible de 
prendre toutes les valeurs possibles entre deux limites extrêmes a et b. Soit p(x) 


2 
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‘ une fonction telle que p(x.)dæ soit le nombre d'expériences ayant révélé une 
valeur de x comprise entre æ, — + dx et x, + + dx, et supposons que dx soit un 
infiniment petit. La fonction p(x) est le taux de fréquence de la variable x. 


Si N est le nombre d'expériences, on obtient, les intégrations se substituant 
aux sommations, 


f p(t)dæ=N. 


Le nombre d'expériences donnant lieu à l'observation de valeurs de x com- 
prises entre a et 6 (a<a<f$<b) est 


[7 (æ) dx. 


Quant à la fréquence relative ou probabilité statistique d'obtenir une valeur 
de æ quelconque comprise entre « et B, elle s'exprime par 


[ p (x) dx 


où 


PU) 


To = 
[ p (œ@) dx 


est le taux de probabilité statistique. 

En pratique, un relevé relatif à N expériences sur le compte d’une variable 
continue se fait dans l'ignorance de la fonction p(x). Pour se rapprocher le mieux 
possible des fréquences théoriques, infiniment petites p(x)dx, on partagera l’in- 
tervalle (ab) en n intervalles égaux, de longueur h, centrés sur les points æ,, æ 
sd sHelSuqUe 


k 

17 RE ln rer 
k 

En = DE 


Dr = Lit (R —1)À. 


On notera les fréquences absolues relatives aux n intervalles ainsi définis par 
o(t)h, o(t2:)h, .…. o(x,)h. 

On définit ainsi une fonction à paliers y=o(x), qui, dans chaque intervalle 
partiel, reste constante et égale à o(æ). Dans ces conditions, la fréquence e(x,)h 
a pour mesure l'aire du rectangle ABCD (fig. 3) et la somme des aires de tous les 


D '/ 
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rectangles analogues, de même base h et de hauteurs p(t), p(r2), .…. ox), est 
égale à N : 
“ n 

D (Cr) À =IN 

21 


y- pl 


FIG. 3. 


$S 2. Problème fondamental de l’ajustement des séries statistiques. 


Si l'on veut tirer des conclusions, susceptibles de s'exprimer en langage 
mathématique, d’un relevé statistique portant sur la mesure d’un caractère jouant 
le rôle de la variable statistique, on ne peut se contenter d’une distribution expé- 
rimentale brute, telle que celle que nous venons de décrire. On ne peut se conten- 
ter davantage de la caractérisation d’une série statistique au moyen de quel- 
ques constantes qu'on estime significatives (moyennes, dispersion, mode, 
médiane, etc.). Ce sont, en effet, des résultats ou bien de l'application de formules 
élémentaires, ou bien de la considération d’un graphique « brut », ou bien 
encore d’une numérotation dans un ordre naturel de croissance ou de décroissance 
des valeurs expérimentales. Ces nombres sont si peu « caractéristiques » qu'il 
_n'est guère compliqué de trouver des séries dont la représentation géométrique 
fait apparaître une complète dissemblance, mais qui peuvent avoir en commun 
certains de ces nombres calculés au moyen de formules simplistes,. 

Si l’on opère des mesures dans des conditions matérielles identiques sur des 
objets, des êtres d'un même groupe ou de même espèce, on constatera que les 
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lois de distribution des fréquences obtenues de plusieurs séries d’un grand nom- 
bre N d'expériences ne sont pas rigoureusement identiques. Les répartitions brutes 
se « ressemblent », elles sont assez voisines les unes des autres, mais on observe 
des écarts d’une répartition à l’autre. 

Il s'ensuit, nécessairement, que les constantes statistiques relatives à chaque 
série de N expériences varieront de l’une à l’autre de ces séries. 

Quel parti peut-on tirer, dès lors, des résultats d’une série de N expériences 
sur des individus pris arbitrairement parmi une collection plus nombreuse ? 
Quelles conclusions l'étude de cette série permet-elle d'établir, conclusions appli- 
cables à toute la collection, à toute l'espèce ? Quels sont les fondements logiques 
de ces conclusions et quelles en sont les limites ? 

Il peut être répondu à ces questions en disant que les résultats expérimentaux 
bruts sont les données d’un problème dont la résolution et les conclusions sont 
fondées sur la théorie des probabilités. Les limites de ces conclusions dépendront 
directement des écarts qu'on pourra considérer avec une quasi-certitude comme 
ne devant jamais être dépassés. 

On recherchera si la série brute étudiée ne présente pas certaines régularités 
qu'on pourrait traduire par une loi mathématique simple. C’est dans la recherche 
de cette loi mathématique que consiste le problème de l’ajustement de la série 
statistique à une variable. Cette loi s’exprimera au moyen d’une fonction de la 
variable statistique. 

Deux problèmes, algébrique et analytique, consistent à déterminer une fonc- 
tion y= f(x) au moyen de la donnée des valeurs prises par cette fonction pour un 
certain nombre de valeurs de +, ou au moyen de la donnée des valeurs de l’inté- 
grale définie de cette fonction sur un certain nombre d’intervalles contigus de la 
variation de x. Ainsi posés, ces problèmes sont complètement indéterminés et il 
convient de choisir à priori la forme des fonctions. 

Or, il ne peut être question d'identifier les valeurs directement observées aux 
valeurs prises par la fonction inconnue, afin de la déterminer par des méthodes 
d'algèbre ou d'analyse mathématique. On sait, par expérience, que les valeurs 
observées fluctuent d’une série statistique à l’autre. 

Ces valeurs observées doivent donc être « compensées » afin de devenir les 
« vraies valeurs », c’est-à-dire les valeurs Les plus probables. En conséquence, on 
tentera toujours d'identifier légitimement, hypothèses et expériences à l'appui, 
la distribution observée à l’une des distributions caractéristiques étudiées de façon 
générale dans la théorie de la statistique mathématique. 

Il ne convient pas ici de se livrer à une étude détaillée des multiples méthodes 
conduisant à une solution convenable du problème. Disons seulement qu'elles 
conduisent toutes à des solutions approchées, solutions dont il est possible de 
pousser, dans bien des cas, l’approximation aussi loin qu'on le désire. 

Notons dès maintenant que, dans le choix des approximations, une juste 
mesure s'impose. Si une approximation grossière n'est pas de mise, une approxi- 
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mation exagérée devient illusoire, car elle revient à apporter des corrections 
analytiques insignifiantes par rapport aux écarts entre les grandeurs issues de la 
mesure directe et les grandeurs théoriques. 


S 3. Moyennes et écarts. Moments. 


Il sera exposé ici, sommairement, les méthodes essentielles susceptibles 
d'application dans cette étude biométrique et dans les études similaires. Toute 
démonstration mathématique a été exclue et il est demandé au lecteur d'accepter 
les résultats énoncés. Néanmoins, il a paru indispensable de définir un certain 
nombre de notions d’usage très fréquent. 


1. Distributions discontinues : 


Utilisant les notations déjà rencontrées au $ 1, nous définirons la moyenne 
arithmétique des x, par 


La quantité x,— A est l’écart absolu de la moyenne arithmétique. 


La dispersion ou écart quadratique moyen & est la racine carrée de la moyénne 
arithmétique des carrés des écarts absolus : 


st — RD (i—AYe (2) 


Le moment d'ordre s de la variable x est, par définition, la moyenne arith- 
métique des s— n% puissances des x, : 


dl n 
NO oO (T0): 
Li D o(æi) 
Le moment d'ordre s de l’écart absolu est, de même, 


RON: (œ; — A} p(æœé). 


DS 
+L 

2. Distributions continues : 
Les mêmes définitions se retrouvent dans le cas des distributions continues, 


mais les signes d'intégration ( sont à substituer aux signes de sommation Y. 


La moyenne arithmétique ou valeur moyenne est 


h 
zh(r) dr [ xT dx 


- 
(72 (2 


et l’on retrouve sans peine l’analogie avec la définition correspondante du 1. En 
effet, p(x)dx est le nombre des valeurs de æ comprises dans un intervalle (infini- 
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ment petit) dx, centré sur le point x; æp(x)dx est la somme de ces valeurs de x, aux 
infiniment petits du second ordre près, et 


[ æp(x) dx 


- 


(14 


est la somme de ces valeurs étendue à tout l'intervalle (ab). 
Enfin, 


À b 
N | x p(x) dx 


est la moyenne de ces N valeurs. 
De la même façon, l'écart quadratique moyen est 


b b 
— [ (æ — A}p (x) dx — [ (x — APT, dæ. 


eo La 
(22 (22 


Les moments M, et x, de la variable x et de l'écart x-A sont respectivement 


1 b 
ME foptam= ( æwr, de 


- 


el dpi? è 
- H=X [ (x — A)'p(x)dx = [ (x — AY T, dæ. 


En particulier, qu'il s'agisse d’une distribution continue ou d’une distribu- 
tion discontinue, on a 
Net MA 
Vo=M=1; mw—=0; jp = 07. 


3. Correction des moments : 


Dans le cas d’une distribution continue, on ne connaît pas, à priori, la forme 
de la fonction T,. Tout au plus peut-on se permettre certaines hypothèses au sujet 
de cette fonction, hypothèses vraisemblables en raison de confirmations très fré- 
quentes. On peut calculer des moments « observés » m, par la formule 


nn 


Ms ai © (%) h, 


N = 


et les hypothèses légitimes sur la fonction T, permettent d'exprimer les moments 
théoriques M, en fonction des moments observés. On est conduit ainsi aux for- 
mules de SxePppArD [R. Risser et C. E. Traynarp (|, d'emploi fréquent. On a 


NE Me 


1 
TR te PO 

il 
M, — Care nu h?, etc. 


(55) RISSER, R. et TRAYNARD, C. E., 1933, pp. 35-37. 
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On ne tiendra compte des termes de correction que pour autant qu'ils ne 
soient pas négligeables vis-à-vis des moments observés. 
4, Signification de la dispersion : 
Dans les deux cas de distributions continues et discontinues, considérons les 
expressions 
b 4 n 
SI= ( (G—X}T,dæeltS — N » (ri — X} 0 (æ:), 
î=1 
&, 


où X est une quantité variable. Les valeurs de X rendant minimum ces expres- 


\ 


: Fe SE . dS : NE 
sions obéissent dans les deux cas à l'équation e —=0, ou bien, en divisant au 
a 


préalable par — 2, à 
| ‘Pb 
G=X)T,dœ=,0 


ou 


ou bien encore 


La moyenne des carrés des écarts à la moyenne arithmétique est plus petite 
que la moyenne des carrés des écarts à tout autre nombre. 

De plus la forme même de ç° indique, et un raisonnement rigoureux du 
calcul des probabilités le démontre, que plus & est petit, plus la répartition est 
concentrée autour de la valeur moyenne, d’où le nom de dispersion. 

Il convient de considérer, non pas tant la valeur absolue de & que sa valeur 
relative par rapport aux x,, ou par rapport au nombre À qui en caractérise l’ordre 
de grandeur. À cet effet, on utilise souvent, comme mesure de la concentration, 
le coefficient de concentration _ On suppose les x, essentiellement positifs, sans 


quoi on pourrait obtenir A—0 et ce cocfficient serait dépourvu de sens. 


5. Moments d'ordre supérieur : 


- Ils jouent un rôle de première importance dans les diverses solutions des 
problèmes posés par l'ajustement des séries statistiques à une variable. Les 
fonctions qui représentent théoriquement les distributions observées sont, en 
effet, très souvent exprimées au moyen de sommes de fonctions simples dont 
les coefficients s’expriment au moyen de ces moments. L'exposé, très succinct, 
du $ 4 en donne une explication. 
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$S 4. Principe des moments. 


Il s’agit d’un principe sur lequel on peut fonder la solution du problème de 
l'ajustement. Il revient à déterminer une fonction T,("), telle que ses n+1 pre- 
miers moments des écarts bo, Li, ..., pe, Soient égaux à ceux des écarts observés 
(éventuellement corrigés au moyen des formules de SHEPPARD). 

On devra donc réaliser les conditions 


L2 


400 +00 
Us = [ aTwda= | LT do S—= 042 27, 
—0 do 


où T, représente le taux théorique, T,(® le taux obtenu par l’application du 
principe des moments. T,°® constitue donc une approximation de T, et, de 
façon générale, cette approximation est d'autant meilleure que le nombre 
n+1 de moments est grand. Cette approximation est, suivant une première 
méthode, obtenue au moyen de développements en série limités à un nombre 
fini de termes. 

Cette première méthode est celle du développement en série de polynômes 
d'Hermire [R. Risser et C. E. Traynarp (), R. Decruegiz (°°)]. Sous certaines 
conditions de continuité et de dérivabilité, toute fonction y=f(f) peut s’écrire 
sous forme d’un développement convergent : 


y = (@o + A Hi + & H+….)e Fr, 


ad" 
Hu EN En (EME) 
est un polynôme de degré n. 
On a 
H,= 1, H=22k'E, H—4k1E2- 292, …. etc. 


Les propriétés analytiques de ces fonctions H, permettent de démontrer que 
le développement dans l'intervalle (—®, +) 


E2 ; 
PERTE emfi-hnt+(h En] (1) 


\'ELTS 


convient à de nombreuses lois de probabilité y=T,. Il suffit de poser £=zx— A 
et d'y remplacer les 4, par les moments des écarts, tirés de l'observation. 

Une seconde méthode, très précise et conduisant à d'excellents résultats, est 
celle de PEARSoN. Elle est basée en premier lieu sur une propriété établie par 


(54) RISSER, R. et TRAYNARD, C. E., 1933, pp. 218-221. 
(55) DELTHEIL, R., 1930, pp. 87-90. 
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PEARSON, revenant à dire que de nombreuses lois de probabilité y=T, obéissent 
à l'équation différentielle 
dy x—a 
— = d æ, 2 
y F(x) @) 


où F(x) est un polynôme de la forme 
F (&) == ao + d4 @ + ae D? + +. 
On a alors, en général, 


DD ENT) 
Der 
ral F (æ) » 


C étant déterminé par le fait que l’intégrale définie de y dans le domaine de x 
est égale à 1 (relation de certitude mathématique). 

On peut démontrer que la connaissance des moments observés permet de 
déterminer &, &, &, &, …, etc. En particulier, si l’on connaît M, et M,, on peut 
déterminer a et «, dans l'équation 


so dæ. (3) 
On trouve, en effet, 
a = M,, 4 — M — M, 
et 
__@—My 
VE Sn 2 (M — ME) (4) 
V27(M,— M) 


C'est la solution du problème dit des deux moments. 

L'utilisation de trois, quatre... moments conduit à des solutions plus géné- 
rales et plus précises. 

On obtient alors des lois ou taux de probabilité dont la représentation gra- 
phique conduit aux courbes de Pearson [R. Risser et C. E. Traynarp (°)]. 


S 5. La loi de GAUSs. 


L'application de l’une des méthodes précédentes à la solution de problèmes 
d'ajustement statistique permet donc d'obtenir des solutions de forme très géné- 
rale, et la précision de la représentation analytique peut être poussée aussi loin 
qu'on le désire, si l’on veut pousser le calcul des moments d'ordre supérieur 
jusqu'à un ordre suffisamment élevé. Notons toutefois que cette précision devient 
tout à fait illusoire à partir de l'instant où la correction apportée par le calcul 


_de moments d'ordres plus élevés est inférieure à l'erreur de mesure et, à fortiori, 


à l'erreur arithmétique commise en négligeant des termes à partir d’un certain 
ordre de grandeur très petit. 


(56) RISSER, R. et TRAYNARD, C. E., 1933, pp. 52-62. 
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D'autre part, on admet couramment, à titre de postulat, que les écarts des 
valeurs mesurées à la « vraie valeur » de la grandeur mesurée sont distribués 
symétriquement autour de cette « vraie valeur ». Autrement dit, on admet des 
probabilités égales pour des écarts de signes contraires, mais égaux en valeur 
absolue. Enfin, si l'on pose que la « vraie valeur » est la moyenne arithmétique 
des valeurs observées, on justifie sans peine la loi de distribution bien connue 
sous le nom de loi de Gauss. | 

Cette loi n’est pas une loi générale de distribution des grandeurs mesurées 
autour d’une valeur remarquable, dénommée « vraie-valeur », aussi remar- 
quable que dépourvue d’une définition précise. Elle est le résultat d’une justifi- 
cation excellente, mais non d’une démonstration, basée sur des hypothèses ayant 
la valeur de postulats, qu'on est libre, en toute logique, de rejeter et de remplacer 
par d’autres, plus complexes. 

L'expression mathématique de la loi de Gauss est la suivante : 


n (x — A} 
Yy=T, = ——2 21540 
ç \2 T 
ou encore 
æ-— A} 
KA on RC 
p(X) = —— 0e 2e 14, (5) 
C3 \2 T 


Les hypothèses sur lesquelles cette loi est fondée, si elles ne manquent pas 
de légitimité, manquent, comme on le voit, de généralité dans bien des cas 
d'application. 

Ce manque de généralité se remarque d’ailleurs aussi dans les exposés som- 
maires des méthodes des polynômes d'HERMITE et de PEARSON, rappelées précé- 
demment. | 

En effet, dans le développement de la formule (1), si l’on néglige les coef- 
ficients de H, et H,, on retrouve 

ju @A} 


VE : € 2 Le ) 
27 fu 


: = MU 
et la solution de l'équation de PEARSON (3), cas particulier de la formule (2), est 
également l'expression de la loi de Gauss. 
On a, en effet, 


lol) 
br) (æ—AÿT,dæ 


b D 
— |  22T; dx — 2 |" x TT; dx + A? 
: à 


= MS AS AE MAO EM MS 


a 
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I convient d'insister fortement sur ce fait : la loi de Gauss n’est qu’un cas 
très particulier des lois de distribution des séries statistiques à une variable. Elle 
na pas le caractère universel de loi des erreurs de mesure qu'une habitude et 
une tradition de l'emploi lui confèrent, d’une façon implicite et quasi incontes- 
table, dans l'esprit de nombreuses personnes. 

Si son emploi est abondamment répandu, c’est que, dans nombre de cas, 
elle donne une bonne représentation de la réalité observée, sans qu'il soit toute- 
fois permis de prétendre que ce soit la meilleure. La loi de Gauss possède quel- 
ques propriétés analytiques qu’il est, en général, très aisé de confronter avec des 
propriétés correspondantes de la distribution observée. Si des écarts trop consi- 


- dérables sont constatés entre la théorie et l'expérience, il faut délibérément 


adopter une autre loi de distribution. 

Dans l'étude de la distribution des indices caractéristiques des populations 
d’'Anthraconauta faisant l’objet de ce travail et exposés au chapitre suivant, il 
a été fait un usage général de la loi de GAUSS, parce que l'ajustement pouvait être 
considéré comme satisfaisant et aussi parce que le nombre réduit des exemplaires 
étudiés (plusieurs centaines) ne permettait pas de tirer des conclusions tendant 
à rejeter cette loi. 

Il nous est apparu cependant que les distributions observées étaient plus 
concentrées que la distribution de Gauss. Nous pensons qu'une étude portant sur 
un très grand nombre de coquilles révélerait une distribution quasi symétrique, 
à concentration plus forte, ce qui s’exprimerait aisément par l'emploi de moments 
jusqu’au 4° ordre dans les méthodes de développement en séries de polynômes 
d'HERMITE et de PEARSON. 


S 6. Propriétés de la loi de Gauss. 


1° De façon générale, nous avons vu que 


LA 
(42 


b (] 
fouwde-x ou [ LÉGLEN, 


l'intervalle (a, b) étant celui dans lequel se trouvent certainement toutes les 
valeurs de x. 
Ici on a 


+0 4 LPS), 
{ 1 Do L = A, 
[ex 


ce qui, géométriquement, signifie que l'aire située sous la « courbe en cloche » 


x 


de Gauss, entre cette courbe et l’axe des x, de —@ à +, est égale à 1. Les 
bornes de l'intervalle où se trouvent certainement toutes les valeurs de x sont 
—o et +oo. Il est bien évident qu’en pratique, les bornes de cet intervalle 
sont finies, mais on remarque qu'à partir d’un certain écart de la moyenne A, 
l'aire située sous la courbe en dehors de cet écart est pratiquement nulle. 
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La probabilité d'obtenir une valeur de æ comprise entre x, et x, est 


DEA 
Pour calculer ces probabilités, on opère la substitution {= V2 
De cette façon, la probabilité précédente s’écrit 
À le 
UE — | exit: ji (6) 


T1 


(7 


On a dressé des tables donnant les valeurs de la fonction de Poisson : 


1 +À 2 À 
O() = — { er ® dt = | eTË dé, 


Vz —À Vz 0 


qui est la probabilité pour que 
A— Xo V2 one A + ho V1. 


Ces tables, établies pour un grand nombre de valeurs de À, permettent 1 
calcul des probabilités (6). 
On constate, notamment, que pour 
1-00) 0e 
—9, A(À) — 0,995, 
, 0 (À) — 0,99998, 
, 0 (X) = 0,99999998. 


(e) 


UN 


On peut admettre que, pratiquement, les écarts de la moyenne arithmétique 


ne dépassent pas en valeur absolue 35 V2, ou bien qu'on a, avec une quasi- 
certitude, 


A3 26 <m<A+3\2. 


2° Définissons à présent l’écart probable et l'écart moyen, éléments carac- 
téristiques de la loi de Gauss. 

L'écart probable E, est celui pour lequel les probabilités d'être dépassé et 
de ne pas être atteint sont égales, donc, toutes deux égales à 1/2. Autrement dit, 
la moitié des valeurs observées, supposées obéir à la loi de Gauss, se trouvent 
dans l'intervalle (A—E,, A+E,). 

Ayant 

pen” ja e-® dt, pour À — 0,4769 
CR Ten 


on à 
E, — 0,4769 V2 & = 0,67456. (D 
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L'écart moyen E,, est la moyenne des valeurs absolues des écarts x— A : 
+0 4 no) 
E, = [ Ga fa 202 dæ:; 
T 


Be 070770. (8) 


Les constantes ©, et E, sont des mesures commodes pour les écarts de la loi 


de GAUSS. 


D'autre part, les relations précédentes entre 5,, E, et E, sont des moyens 
de contrôle du bien-fondé de l'identification d’une répartition expérimentale à 
cette loi de distribution. 


$S 7. De la moyenne arithmétique considérée comme 
une variable statistique. 


Soit une série de N expériences portant sur une variable statistique x et 
ayant fourni les valeurs x,, æ, ..., æ,, certaines valeurs pouvant se reproduire 
plusieurs fois. 


n? 


Posons z= T 2 * ; chaque série de N épreuves sur le compte de x 
fournit des valeurs x,, x,, ..., æ, et une certaine valeur de z, z est une varïable 
statistique au même titre que les æ,, æ,, ..…, æ,, et nous l’appellerons la variable 
moyenne. 


Affectant des indices x et z les symboles des moyennes et dispersions, il est 
aisé de démontrer que 


gs? — LE 
3 N°’ 


ou bien 


Sœ 
Ty —= ——:. 
Vx 
Si l’on connaît le taux T,, on peut calculer T, et, notamment, on peut 
démontrer que si la variable x obéit à la loi de Gauss, il en est de même de la 
variable moyenne z. On obtient alors, de façon précise, 
(is AP 


(10) 


ns us 1) 


Si æ ne suit pas rigoureusement cette loi de distribution, T, n'y obéit pas, 
mais la formule [10] montre que si N croît indéfiniment, 6, a pour limite 0. 
Ceci suffit pour prouver que si le nombre d'épreuves croît indéfiniment, la 
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variable moyenne z est « concentrée » en A(x), ce qui signifie que l’erisemble 
des valeurs quasi certaines de z est constitué par un intervalle tendant vers 0 et 
centré sur A(æx). 

Dans ces conditions, que nous qualifierons d’asymptotiques, il est très 
approché d'admettre que z obéit à la loi de Gauss [11]. 
| Dans le cas de la présente étude, cette approximation est excellente et 
d'autant meilleure que N est plus grand. 


On peut alors admettre que l'intervalle des valeurs a certaines de z est 
défini par 


AGE INSEE nt do (42) 
VX 


AVES 
3 V2 ere | (15) 


constitue, en valeur absolue, une limite supérieure de l'erreur commise sur la 
moyenne. | 


et que 


Cette conclusion est d'une grande importance pratique, comme nous aurons 
l'occasion de le voir. 


Comme au $ 6, nous obtiendrons les expressions de l'écart probable sur la 
moyenne E,°” : 


[°3 ’ 
ES) — 0,6745 —Z., (14) 
. VW 
et de l’écart moyen sur la moyenne 
EP = 0,7977 <=. (45) 
N 


8 8. Symétrie et asymétrie d’une distribution observée. 


La notion intuitive de la symétrie d’une distribution est très claire : une 


distribution est symétrique si les valeurs de la variable étudiée se répartissent 
symétriquement autour d’une valeur particulière, 


Analytiquement, la distribution de la variable x est symétrique autour de a 
si la fonction T, est paire en x — a : 
on f[(æ KE a} |, 
c'est-à-dire si T, ne s'exprime qu'au moyen des puissances paires de x— a, et 
réciproquement. On démontre alors que A=a. 
Géométriquement, la courbe y=T, est symétrique par rapport à l’axe x= A 


(fig. 4). 


Une propriété remarquable s'en dégage : tous les moments des écarts d’ordre 
impair d’une distribution symétrique sont nuls : 
q 


== Us =. long = 0. 


LLR. 
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La loi de Gauss, en particulier, est une distribution symétrique. 11 est aisé 
d'obtenir une forme analytique générale de toutes les lois de distributions symé- 
triques dans l'intervalle de + à —. Il suffit, dans le développement en 
série de polynômes d'HERMITE, d'annuler tous les termes contenant lu, (ts, 
Hess +. Monyie , 

En pratique, il est rare, parce que très peu probable, d'obtenir une répar- 
tition observée rigoureusement symétrique. En effet, même si l’on sait, à priori, 
que la répartition est théoriquement symétrique, ou si l’on a des raisons légi- 


y=T,- [ æ-AŸ] 


A 


FIG. 4. 


times de le supposer, la loi des écarts dus au hasard ou aux erreurs joue toujours 
lorsqu'on tente une réalisation expérimentale. On ne peut donc affirmer l’asy- 
métrie d’une distribution sur la foi de l'observation montrant une asymétrie 


modérée. En fait, l'hypothèse légitime décide de la symétrie et l'hypothèse ne 


peut être infirmée que par des divergences flagrantes. 

La meilleure mesure pratique est la constatation de la non-nullité des 
moments Le,s1, OU, du moins, du fait que ceux-ci sont, en valeur absolue, d'ordres 
de grandeurs comparables aux moments d'ordre pair. 

Les mesures empiriques de l’asymétrie (« skewness » des Anglo-Saxons) ont 
un caractère synthétique; mais, comme d’autres constantes statistiques, elles ne 
donnent qu'une indication générale imprécise. On peut trouver un exposé 
détaillé de la méthode de ces mesures dans le traité remarquable de G. Unxy 
YuLe et M. G. Kenpazz (°, 


(57) Upny YULE, G. et KENDALL, M. G., 1945, pp. 96, 98, 162-164. 


30 ÉTUDE BIOMÉTRIQUE DES ANTHRACONAUTA 


8 9. Méthode pratique de calcul des moyennes et de la dispersion. 


Soit une variable statistique continue æ susceptible de prendre toutes les 
valeurs d’un certain intervalle, Supposons que celui-ci soit, pour les nécessités 
du relevé des fréquences, partagé en intervalles partiels de longueur h. 

Comment fixer le choix de h ? Si l’on désire se rapprocher des exigences 
théoriques, il faudrait que h fût infiniment petit. D'autre part, il faut aussi que 
l'intervalle considéré contienne des valeurs de æ donnant lieu à des observations, 
car le nombre des valeurs observées de x est forcément fini. On prendra donc 
un nombre d'intervalles dépendant de l'intervalle total, ni trop petit, pour que 
la courbe à paliers garde une signification, ni trop grand, pour qu'on puisse 
effectivement conserver des valeurs de # et pour donner plus de régularité à ce 
qu'on observe. Il faut noter que ces intervalles doivent, au moins, être d’une 
grandeur double de l'erreur de mesure. 

En général on prendra de dix à vingt intervalles. 

Les calculs de moyennes, de dispersions et, en général, de moments peuvent 
être très longs, si l'on emploie simplement les formules du $ 3. Aussi opère-t-on 
souvent une substitution de variable très simple, ramenant à l'unité l’intervalle 
de longueur h. On passe de la variable x à une variable X par la relation , 

DA ENT, (16) 
où +, est une valeur de x convenablement choisie, voisine, à première vue, de la 
moyenne arithmétique, de sorte que X prendra des valeurs positives et négatives, 
assez petites, dont la moyenne sera voisine de 0. 

On voit que si &, — x, = h, on à R(X,;, — X,) = het X, — X, = 1. 

L'intervalle h est donc bien ramené à l'unité. 


Cela étant, désignant par A(x), 6,, A(X), 54 les moyennes et écarts quadra- 
tiques de x et X, on a immédiatement 


A(x)=hA(X) +: (17) 


Soustrayant [17] de [16], on à &æ— A(x)—h[X — A(X)] et 6? — h? 5%,ce qui 
ne dépend pas de «,. 

Ces formules sont d’un usage facile et nous en ferons un emploi fréquent. 

Par exemple, soit une variable statistique x observée dans l'intervalle 
118,5 < x < 157,5, soit h — 4, d’où le nombre d'’intervalles partiels n = 11. 


On a 2,=118,5+ À —115,5, 2—119,5, #,—198,5; … 210 151,5, œ — 155,5. 


Posons æ = 4X + 135,5. On trouve alors 


Me X; = —1 X, = 3 
X, = —4 Xe —0 X,9 = 4 
X3 = — 32 NOTE XX —5 


X, — —2 Xe 
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Soit, enfin, le relevé des fréquences : 


SCD —2 . eo (ms) = 37 (a) h = 5 
p (Xe) h="2 © (%e) h = 45 @ (io) R — 
o(æ)h =7 © (æ@r) h = 27 o (æu) h = 1 
p (2) h = 14 o (æ3) h = 18 


al 


D œ(a) h=—159. 


On a 


A(x) — NES ER CDS) 4x (—2) 87 x (1) 


17 


1 
LH 4BX 0 L27X 1 HISX2HLBXFTIXATIX que = 0,106 
a) 
et 


A (D) = 4X (—0,1) + 195,5 — 135,1. 


: 


Calculer les carrés X, — A(X) est fastidieux. On dispose d'une autre relation 
simplifiant encore le calcul de la dispersion. Ayant en général 


Fb 
s,— | (x — APT, dx, 


a 


ou bien 
bd 
Lac | (@—2A% + ADT, dæ, 
a 
il vient 
, ro Pb bd 
Ge | da T, dx —2A ] æT,d4x+.A2 | T,.dæ. 
(2 (22 a 
6% — M, — 2 A? + À?, 
6% = M,— A? — A(x?) —[A(æ)P. (19) 


En particulier, 


4 
GE = A(X) —[A(X)P = [2 X 25 +2 x 16+7X 9414 X 4 


+ 37X1+27X1+18X 4+5X 9 +1X 16 + 1X 25]— 0,0112 — 2,6484 
et 
Cx 11,02 
4 102048, 


- S 10. Liaisons entre deux variables statistiques. Corrélation. 


Supposons que, dans un ensemble d'objets de même espèce, la mesure d’un 
caractère bien défini soit désignée par x, et celle d’un autre caractère par y, ces 
caractères étant supposés linéaires. Pour chaque spécimen on obtiendra ainsi un 
couple de grandeurs (x, y) et la question à résoudre, exprimée sous une forme 


3 
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générale et toutefois imprécise, est la suivante : existe-t-il une liaison entre les 
caractères x et y mesurés sur chaque objet ? 

Pour y répondre, certaines notions d'analyse mathématique et de calcul des 
probabilités doivent être rappelées et précisées. 

Considérons, de façon générale, une collection d'objets ou d'individus d'une 
espèce, sur lesquels nous mesurons deux caractères, x et y. 

Il peut se faire qu’une valeur x fixe y de manière univoque et réciproque- 
ment. On dit alors qu'il y a une liaison ou dépendance fonctionnelle entre x et y 
comme entre l'intensité du courant et la tension dans un circuit électrique de 
résistance invariable. | 

Il peut se faire aussi que cette dépendance soit plus ou moins lâche, en ce 
sens que pour chaque valeur de æ ou chaque petit intervalle de valeurs de x on 
observe une répartition de valeurs de y ne permettant pas de conclure à une 
liaison fonctionnelle. Précisons cette notion par des exemples bien connus : Si 
l’on mesure la longueur et la largeur d'œufs de poule d’une espèce bien détermi- 
née, le sous-croupe qu'on obtient en fixant la longueur, ou en considérant un très 
petit intervalle de longueurs, n’est pas composé seulement d'œufs ayant une lar- 
seur bien déterminée, ou comprise entre des limites très rapprochées, mais il y a 
aussi, pour chaque longueur, une répartition des largeurs suivant une loi de fré- 
quence rappelant la loi de Gauss. Il en est de même si l’on mesure la taille et la 
largeur d’épaules d’un groupe d'hommes assez nombreux. Le sous-groupe des 
hommes de même taille n’est pas composé d'hommes ayant, de surcroît, la même 
largeur d’épaules; mais on constate également une répartition d’allure symé- 
trique rappelant la répartition suivant la loi de Gauss. 

On dit alors qu'il y a liaison stochastique entre les deux variables statistiques, 
longueur et largeur des œufs ou taille et largeur d’épaules, et d’une façon géné- 
rale’eritre æet.y; 


Mais deux cas encore sont à distinguer : 1° la répartition des largeurs des 
œufs (ou des largeurs d'épaules) ne change pas, quelle que soit la longueur des 
œufs (ou les tailles des hommes); on dit alors que les deux variables x et y sont 
indépendantes au sens des probabilités; 2° la répartition des largeurs dépend de 
la longueur, et dans ce cas x et y, longueurs et largeurs, sont en corrélation. 


Ces notions ne peuvent être exprimées avec précision que sous la forme ana- 
lytique. 
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$S 11. Distribution des fréquences d’un complexe 
de deux variables statistiques. 


Il convient, pour avoir une idée plus exacte de la corrélation, de poser quel- 
ques définitions et de détailler l’origine des nombres par lesquels on tente d’en 
donner une mesure. 

Considérons le cas de deux variables, x et y, continues, analogues aux lon- 
gueurs et largeurs d'œufs ou aux tailles et largeurs d’épaules, en négligeant les 
cas où les deux variables sont discontinues, ou bien, où l’une est continue, l’autre 
discontinue. À chaque mesure, sur un même individu, des deux caractères x et y, 
faisons correspondre un point d’un plan rapporté aux axes O,, O,. Nous obtenons 


ÿ 


F1G. 5. — Liaison fonctionnelle. 


ainsi un nuage de points dans le plan. Si la liaison est fonctionnelle, tous les 
points sont situés sur une même courbe ou au voisinage immédiat de celle-ci. 
Dans la négative, on peut observer une répartition des points avec une concen- 
tration plus ou moins prononcée autour d’une ligne. Si cette concentration est 
très prononcée, la liaison stochastique est proche de la liaison fonctionnelle 
(fig. 5 et fig. 6). Autrement, en découpant, dans le plan, une tranche de largeur 
suffisamment mince, perpendiculaire à O,, on observe une certaine répartition 
de fréquence de y, variable en général, avec l’abscisse de la tranche choisie. 
C'est le cas de la figure 6. Si cette répartition ne varie pas, æ et y sont indépen- 
dantes (fig. 7). 

Traduisons analytiquement ces notions intuitives. 

Supposons que les couples de valeurs observées (x, y) définissent des points 
situés dans un domaine D du plan xOy. Couvrons ce plan d’un quadrillage, et 
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soient h et k les côtés des rectangles parallèles aux axes O, et O,. Soient +,, &, les 
coordonnées du centre d’un quelconque des rectangles obtenus (fig. 8). 
Désignons par Ÿ(x,, y)h.k le nombre des points du rectangle de centre (x,, y,) 


ÿ 


x 
x 
x] * 
x] * 
x[x x 
x a + La 
x X* 
AR XXe k 
TS x X 
Le x 
RAR IARATS LL x 
EAUX XX * 
IXXXEXS X x x LS 
XX XX x ST 
x 
y. 4 
x 
O - 2e 
h 
F1G. 6. — Corrélation. 
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F1G. 7. — Indépendance au sens des probabilités. 


effectivement dénombrés au cours du relevé. Si N individus ont fait l’objet de la 


mesure des deux caractères æ et y, on a évidemment, les sommes X, et à, étant 
étendues à toutes les valeurs de i et j, 


SE, Lx, y5)h, k = N. 
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On définit ainsi une fonction discontinue z=4{x, y), dont la représentation 
dans un système de trois axes rectangulaires est une surface constituée par la 
juxtaposition de parallélipèdes droits, de hauteurs d(x,, y,) et de bases égales, de 
côtés h, k (fig. 9). 


FIG. 9. 


Ici également se pose, de façon générale, le problème de l'ajustement des 
distributions à deux variables. Il peut être résolu par des méthodes issues du 
calcul des probabilités, et sa solution revient à définir une fonction p(x, y) ou une 


fonction T,, FAR en principe continue, telles que 


un 
[fre y)dæ dy = N, ou Î| T, dæ dy = 1, 
+ D 


D 
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et que 


[lr (æ, y) dædyret [| THdvay 
C8 NS eu F 


soient les nombres des points observés dans le domaine F, compris dans la limite 
de D, et la probabilité statistique d'obtenir un point compris dans les limites du 
domaine F (fig. 8). 

La fonction T,, est le taux de probabilité des deux variables x et y, dans le 
domaine D. 


8 12. Dépendance et indépendance au sens des probabilités. 


Nous avons déjà vu, au $ 1 et au $ 2, de quelle façon peut se présenter la 
liaison entre les caractères æ et y, et nous avons utilisé les termes « dépendance 
fonctionnelle » pour qualifier la liaison entre x et y, s'exprimant par une rela- 
tion Y=/1@). 

Il s’agit d’une dépendance s'exprimant sous la forme d’une fonction au sens 
de l’analyse mathématique. 

La dépendance au sens de la théorie des probabilités s'exprime d’une façon 
différente [P. Van DEUREN (°*)|. 

Considérons deux événements, À et B, el un événement constitué par 
l’arrivée simultanée des deux événements À et B. Désignons par AB cet événe: 
ment « composé ». Un principe fondamental du calcul des probabilités s'exprime 
par la relation 

Pas = Pa PB:a 
— PB Pa:p (19) 
où pag est la probabilité d'arrivée de AB; p,, ps étant les probabilités d'arrivée 
de À et de B; p\.8 la probabilité d'arrivée de À si B s’est produit, et pg.a la 
probabilité d'arrivée de B si À s’est produit. 

Si À représente la constatation d’une valeur de + comprise dans l’inter- 
valle (x, æ+dx), et B la constatation d’une valeur de y comprise dans l’inter- 
valle (y, y+ dy), AB sera la constatation d’un couple (x, y) définissant un point 
intérieur au rectangle de côtés dx et dy, et l’on aura 


Le =: Le © nee 
D. (20) 
Les taux de probabilités T,.. et T,. sont les taux liés de y en x et de x en y. 
T,.. dy est la probabilité d'obtenir une valeur de y entre y et y+dy, 
sachant que x est situé dans l'intervalle de x à x+dæ. Autrement dit, T,., est la 


loi de répartition des valeurs de y correspondant à la valeur x quelconque entre 
z et t+ dx. 


(55) VAN DEUREN, P., 1934, pp. 20, 73, 79-81. 
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En général, T,., dépend de x et de y, c’est-à-dire que la loi de répartition 
de y en x varie avec x. On dit alors que x et y sont dépendantes au sens des 
probabilités. 

Si, quelle que soit x, la loï de répartition de y en x reste invariable, autre- 
ment dit, si 

RE: (21) 


æ et y sont indépendantes au sens des probabilités. 


(4 


ETG- 10: 


Notons que la relation [21], en vertu de [20], entraîne 
DUGENTIE 


La dépendance fonctionnelle, analytique, et la corrélation entrent toutes 
deux dans la définition de la dépendance au sens de la théorie des probabilités. 


S 13. Moyennes liées. Lignes de régression. 


Une valeur de x étant fixée, l’ensemble des valeurs de y, observées en même 
temps que cette valeur de x, est contenu dans un certain intervalle D,,, dont 
les extrémités dépendent de x, comme le montre intuitivement la figure 10, et 


l’on a 
(ne dy À, 


Dy:x 


relation de certitude, exprimant que toutes les valeurs de y, x étant fixée, sont 
dans l'intervalle D... | 
Les moments de la loi de répartition T,.,, sont les moments liés de y en x 
et, parmi eux, $ 
ARENA. 2 dy (22) 
Dy:x 


est la moyenne liée de y en «. 
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À chaque valeur de æ correspond ainsi une moyenne liée de y en x, et le 
lieu géométrique des points de coordonnées (x, À,.,) est la ligne de régression 
de yen 

Cette ligne est telle que la somme des carrés (y— A,.,)*, pour toutes les 
valeurs de x et y, est minimum. 

On définirait de même la moyenne liée et la ligne de régression de x en y 
par 

RE Pen e (23) 


Dæ:y, 


Les expressions [22]-et [23], qui sont des fonctions de x et y respectivement, 
sont les équations des lignes de régression de y en x et de x en y. 


S 14. Ecarts quadratique et rectangulaire. Coefficient de corrélation. 


Si l’on ne considère que la seule variable x, indépendamment des valeurs 
correspondantes de y, on peut étudier la loi de distribution de x, ses moyennes 
et moments, et l’on démontre notamment que 


Ca 


À = [| D Toy dx dy (24) 
D 
el, de même, que l'écart quadratique de x est défini par 
= ( AT. de dy. 
; - D 


En fait, si l'on se réfère aux définitions données au $ 3, on constate, sans 
entrer dans la voie des démonstrations analytiques, que AÀ,, calculée selon [24], 
est bien la moyenne de toutes les valeurs de x, la double intégration étant étendue 
à toutes les valeurs de x du domaine D. 

En pratique, nous aurons à calculer À, et 6,, alors que nous sommes dans 
l'ignorance, à priori, de la forme de T,,. Substituant à T,, la fonction discon- 
tinue d(+,, y,) définie par les valeurs (observées) qu’elle prend dans les rectangles 
de centres (x, y,) et de côtés h et k, nous aurons 


4 
A N 2, 2, x; (y (dG:,Y5)h.k, (25) 
à : 
o, — . LE, Z,(a; — A) d(a:, y h.k. (26) 
En posant 
p (25) = E,% (ai, 7) k, 
on a 
Aer ANS: h 27 
Ag = it 0 (2) 9 ( ) 
1 
5? — N », (@; ET A) e (;) h, (28) 
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… g(x)h est le nombre de valeurs de x observées dans l'intervalle CCE 5) 
quelle que soit y. 
De la même façon, 


a = fur de ay, (29) 
D 


= | [ (yre À) T.. d& dy. (30) 
D 


- 


et, en pratique, 


di 
À, — RAY y)h.R. 
d La a 
K R 2197 (5) le (54) 
it 
2 1 DS) 2 fi 
al 


= Nu — Ar (y) k, (82) 
en posant r(y,) = 2x, y)h. 
Outre les écarts quadratiques 6”, et 5°, s’introduit, par la présence des deux 


variables x et y, un « moment » d’une autre espèce, qu’on a appelé l'écart 
rectangulaire : 


ES { (A0) AT, de dy, (83) 
- D 
que nous calculerons par la relation 
AS a 
E— N 2, (@;— A) (y: — A) (2, y) h.k. (34) 


L'exposant rappelle que E?,, est une grandeur de même nature ou de même 
« dimension », au sens de la physique, que 6”, ou 5°, E”,, est donc la valeur 
moyenne du produit (æ—— A.) (y— A,). Il peut être positif, nul ou négatif et, en 
valeur absolue, 
LAC 


ou bien 
Ey 
a | 
Per 
; E? à RTE 
: Autrement dit, le rapport p,,= #%__ [35] peut prendre, à priori, toute 
EC 


valeur comprise entre —1 et +1. Ce rapport est le coefficient de corrélation des 
variables statistiques x et y. 

Supposons que 5, et o, ne soient pas nuls. Dans ce cas, p,,—0 revient à 
Exp) On peut démontrer alors que, si les variables x et y sont indépendantes, 
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Ty 
toujours vraie. Elle l’est dans un nombre élevé de cas particuliers, à vrai dire 


les plus fréquents. 

De plus, on démontre que si |o,|=1, il existe une relation fonctionnelle, 
linéaire, entre x et y du type ax+by+c=0, et réciproquement. 

Plus exactement, cette relation s'écrit, dans ce cas, 


E?. et le coefficient de corrélation de xy sont nuls; mais la réciproque n'est pas 


Y— y = e (x — As). (36) 


æ 


Pour toutes les valeurs de 9, telles que 
A RQ PA RE 


les variables x et y sont dépendantes au sens des probabilités. 


8 15. Signification du coefficient de corrélation. 


Les conclusions qu'ont peut tirer du calcul de 5,,, au moyen de la formule 
de définition [35] et des formules [28], [32] et [34], ne sont précises que dans un 
cas : celui où lo, =1 ou a une valeur si voisine de 1 qu'on peut estimer, en 
raison des écarts possibles, qu’il en est ainsi. Dans ce cas, x et y sont liés par 
une relation linéaire. 


Si p,,—0, ou si 0 <}p,,| 1, on ne peut tirer aucune conclusion précise. 

Parler de corrélation positive lorsque 0 < 0, 1, et de corrélation négative 
lorsque —1<o,,<0, est affaire de définition; parler de corrélation forte, 
lorsque »,, est plus près de 1 que de 0, et de corrélation faible dans le cas inverse, 
est d’une imprécision totale. 

Encore une fois, le problème général qui se pose ici est celui de l’ajuste- 
ment, plus complexe, de lois de répartition à deux variables statistiques. 


8 16. Une loi particulière d'ajustement de répartition 
à deux variables : la loi de BRAVAIS. 


On a ainsi exposé d’une façon sommaire et générale les bases de l’étude des 
distributions de couples de variables statistiques. Il en ressort que toute loi de 
probabilité à deux variables statistiques, bien définie par la donnée du taux T,, 
et du domaine certain d'observation D, donne lieu à des valeurs précises des 
écarts quadratiques et rectangulaires, donc à une valeur du coefficient de corré- 
lation 0... 

Inversement, le coefficient de corrélation ne caractérise pas une distribution 
de couples (x, y), à moins que, pratiquement, on ne trouve que 0,, est nul ou 
égal à un, en valeur absolue. Si l’on a une tendance à donner au coefficient de 
corrélation une signification plus large que celle qui lui est accordée ici, c’est 
que, dans certains cas particuliers, la loi de distribution peut être représentée au 
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moyen d'une fonction à deux variables, où le coefficient 9, apparaît explici- 
tement. 

La loi, généralement représentée par des fonctions de ce type, est la loi de 
BRAVAIS : 


4 1 (he Lx ie = sd 2 Py (Az) (y-A,) Se 
Ne ope | ; k 
La = 2r 6» 0, Vie : € ?2(1-pa) So Se" Ty 5 (37) 
ay 


Entre la loi de Bravais et les lois plus générales fournies par des méthodes 
d'ajustement existent des liens analogues à ceux existant entre la loi de Gauss 
et les lois générales des distributions de la fréquence d’une variable statistique, 
obtenues par les méthodes de PEARSON ou des séries d'HERMITE. 

Cette loi s'applique avec succès dans les cas où les mesures portent sur des 
couples de variables x, y dont chacune suit une loi de répartition rappelant la 
loi de Gauss. Il en serait ainsi pour un des exemples dont nous avons parlé : 
celui de la mesure des longueur et largeur d'œufs de poule d’une espèce. Mais 
dans l'étude qui nous occupe, la loi de Bravais est inapplicable, car les carac- 
tères æ et y sont répartis, séparément, de façon tout à fait quelconque, puisque 
ce sont des mesures de longueurs sur des individus pris au hasard, à tous les 
stades de leur développement. 

Il est aisé de démontrer que les lignes de régression de cette loi sont des 
droites. La droiïte de régression de x en y a pour équation 


0, 
ET à, Ne) " (æ CT A%), (38) 
et la droite de régression de x en y, 


B— A Pay — EG— A) (39) 


8 17. Droites de régression d’une distribution quelconque. 


Les lignes de régression de y en x et de x en y : 


VE Aya 
DIU, 


sont donc, en général, distinctes des droites d'équations [38] et [39]. 

Néanmoins, ces droites, quelle que soit la loi de distribution des couples 
(x, y), conservent une signification. La droite de régression de y en x est une 
droite telle que la somme des carrés des différences entre les ordonnées observées 
de y et celles correspondantes sur la droite soit minimum. La propriété est ana- 
logue pour la droite de régression de x en y. 

Les équations de ces droites peuvent être obtenues en utilisant une méthode 
générale, dite des moindres carrés, sur laquelle on ne peut s'étendre ici [cir. 
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R. Decraæis (*°)]. Cette méthode permet de trouver les équations [Yuze, G. U,, 
et KenpAzz, M. (‘°)], de forme analytique quelconque, de toutes les courbes ayant 
cette propriété. Parmi ces courbes figurent notamment les lignes de régression 
elles-mêmes. 

De plus, si, pratiquement o,,=1, les deux équations [38] et [39] sont identi- 
ques et expriment la dépendance linéaire fonctionnelle qui existe entre x et y: 


hu . (X— As). (40) 


æ 


8 18. Erreur systématique et écart quadratique du coefficient 
de corrélation. 


En toute rigueur, on a 


A 
Pay rest 
ou bien 


ff aa 7. dr dy 


Pay — 


2 


fÉe—a, dedy.[f 7. dx dy 
D D 


T,, étant le taux de probabilité de x et y. 

En pratique, le coefficient de corrélation est calculé empiriquement, en 
utilisant les résultats bruts d’une expérience portant sur N couples de valeurs 
des variables (x, y). Aux écarts rectangulaires et quadratiques théoriques, on 
substitue des écarts qui seront appelés empiriques. On peut établir des formules 
de correction des écarts, donnant la valeur des écarts théoriques en fonction 
des écarts empiriques. Les différences entre écarts théoriques et empiriques 
constituent des erreurs systématiques sur ces écarts, au sens habituel qu'on 
accorde à cette expression, dans la théorie des erreurs, Il en résulte une erreur 
systématique sur le coefficient de corrélation qui s'exprime sous une forme 
générale rendant les calculs extrêmement longs. Cette forme se simplifie consi- 
dérablement si la distribution des couples (æ, y) suit une loi « normale » de 
Bravais. On peut alors considérer que l'erreur systématique du coefficient de 
corrélation est l'erreur par défaut 


ee. 


qu'il convient d'ajouter au coefficient empirique de corrélation pour obtenir de 
façon très approchée le coefficient théorique [R. Risser et C. E. Traynarp (°)]. 


(5°) DELTHEIL, R., 1930, pp. 99-103, 115-134. 
(5°) YULE, G. U. et KENDALL, M., 1945, pp. 209-210, 311-313, 
(5) RISSER, R, et TRAYNARD, C. E., 1933, pp. 195-197, 


_ sois ii 
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De même, l'écart quadratique du coefficient de corrélation, qui est la valeur 
moyenne du carré des écarts entre les coefficients de corrélation empirique 
(considéré comme une variable statistique) et la valeur moyenne du coéfficient 
de corrélation empirique s'exprime également sous une forme générale assez 
rebutante. 

En première approximation, dans le cas d’une distribution suivant une loi 
de BRAVAIS, on a 


> _(A—Ppa) 
Pay N ; 
ou bien 
À — Pay 


Si la distribution des couples (x, y) n’obéit pas à la loi de Bravais, cette 
dernière expression fixe avec une approximation très satisfaisante l’ordre de 
grandeur de l'écart quadratique. 

Ce qu'il faut en retenir, d'une manière générale, c'est que l'erreur systéma- 
- : El 2 : SE El è 
tique est proportionnelle à — et l'écart quadratique à -——-. Pour autant que N 

f N 
soit suffisamment grand, l'erreur systématique devient négligeable devant l'écart 
quadratique. 

Par exemple, supposons 0,,—0,9. On trouve successivement : 


a) Pour N—150 : 


5, — — 0,055, 
12,24 
0,19 
ne 0,9 — 0,0006: 
300 fe 
b) Pour N—9295 : 
0,19 
5, — —— — 0,0127, 
1 
p— %1 \ 0,9 — 0,0004; 
450 
€) Pour N—400 : 
19 
= DE — 0,0095, 
20 
E — NE X 0,9 = 0,0002. 
800 


On voit que l'erreur systématique affecte tout au plus la troisième décimale 
du coefficient de corrélation, puisqu'elle est inférieure au 0,001. Encore s'agit-il 
ici d'expériences portant sur des nombres assez réduits d'individus. Comme 
l'erreur systématique du coefficient de corrélation est une erreur par défaut, on 
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peut en tenir compte en conservant trois décimales du coefficient de corrélation 
empirique, la dernière étant forcée. 

Par contre, l'écart quadratique du coefficient de corrélation affecte très 
souvent la deuxième décimale et l’on peut admettre que le coefficient théorique 
soit situé dans un intervalle dont l’ordre de grandeur est 0,03, centré sur la 
valeur empirique obtenue au moyen de nos procédés de calcul. Get ordre de 


grandeur n’a rien d’absolu; il convient au problème particulier que nous avons 
à traiter. 


C. — MENSURATIONS DES INDIVIDUS. 


Parmi les gisements ayant livré des Anthraconauta belgica cités dans la liste 
aux pages 11 et 12, deux d’entre eux ont fourni un matériel suffisamment abon- 
dant et en assez bon état de conservation pour se prêter à des mensurations (voir 
planches A et B. Ce sont : PI. Chenée, Angleur (5), PI. Andenne, Bois de Stud; 
ils seront désignés par la suite : À. — Angleur et À. — Stud. Pour la compa- 
raison des espèces provenant de ces gisements et attribués à A. belgica, avec des 
espèces d’A. minima, il a été examiné des matériaux de cette dernière prove- 
nant du charbonnage de Wérister, siège José, étage de 242 m., Bacnure Nord, 
stampe des veines Venta à Beaujardin, première veinette sous Grosse, niveau 
n° 48. Ils seront désignés par la suite : À. — Wérister. 

Les mesures efféctuées suivant les méthodes indiquées pages 13 et 14 ont 
porté sur un nombre d'individus variant de 421 à 433 pour A. — Wérister, 193 


à 209 pour A. — Angleur, 146 à 150 pour A. — Stud, suivant les caractères 
mesurés. 


CHAPITRE IL. 


D. — ETUDE DE LA CORRELATION DES CARACTERES. 


INTRODUCTION. — Il a été vu au chapitre IT que la notion de corrélation 
permet d'étudier la variabilité des caractères mesurés, en fonction du dévelop- 
pement (ou âge) de ces individus. Il a été exposé la signification de la corré- 
lation et du coefficient par lequel on peut en donner une mesure, imparfaite 
d’ailleurs, sauf dans les cas extrêmes où il est nul ou égal à +1. Dans ce dernier 
cas, on peut conclure, en toute rigueur, à une liaison fonctionnelle linéaire entre 
les caractères rapprochés. 

Il est évident — et l'attention a déjà été attirée sur ce sujet — qu'une 
distribution expérimentale de couples de variables statistiques, dont on sait à 
l’avance qu'elles sont liées par une relation linéaire, ne donnera un coefficient 
égal à +1 que de façon tout à fait exceptionnelle. En statistique, l’écart est la 
règle et non l'exception. Dès lors, compte tenu des écarts, obtenir un coefficient 
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dont la valeur absolue est de l’ordre de 0,90 est un résultat qui permet de con: 
clure à une dépendance quasi linéaire, surtout si l’on dispose de collections rela- 
tivement peu nombreuses. 


Méruopes. — Les calculs qui suivent ont été effectués selon les principes et 
méthodes énoncés précédemment. 

Afin de calculer les coefficients de corrélation, des tables de corrélation à 
double entrée ont été établies. Pour arriver à ce but, il a fallu : 


1° Partager en intervalles partiels de longueurs convenables les intervalles 
de variation des deux variables mises chaque fois en cause. 


2° Compter les fréquences d'apparition des couples (x, y) dans chacun des 
P 0 
rectangles ainsi obtenus, fréquences désignées au chapitre IT par la notation 
8 ; 8 P E 
o(x,, Y,)hk. On obtient ainsi un relevé statistique brut de la distribution (quel- 
conque) de fréquence des couples (x, y). 


3° Calculer les moyennes liées de l’un des caractères par rapport à l’autre. 


Pour préparer le calcul définitif des coefficients de corrélation et de régres- 
sion, il a été calculé successivement dans chaque cas : 


a) les valeurs moyennes des caractères mesurés; 
b) les écarts quadratiques ou la dispersion des caractères; 


c) les écarts rectangulaires, valeurs moyennes des produits des écarts à la 
moyenne. 


Voici le détail de toute la suite des opérations conduisant à l'établissement 
de ces valeurs : 


Pour chaque population, la corrélation entre les caractères d et a, d et b, 
d et g a été étudiée, le caractère d étant choisi comme élément de référence en 
tant que mesure du diamètre de croissance maximum de la coquille. 

On a désigné par d, ou A(d), & ou A(a), b, ou A(b), g, ou A(g) les moyennes 
des caractères mesurés; par 0,4, 5, 6, 0, leurs écarts quadratiques; par E°,, 
E°», E°,, les écarts rectangulaires. Les coefficients de corrélation se définissent 
alors par 


° 
_ E& 
Paa Ta reg 
Sa * Sa 
2 
: Ed 
GA , 
Sa° 5 
2 
ee 
(2 1 SCENE 
Ta . Sy 


Il a paru nécessaire de donner quelques indications au sujet des nombres 
figurant dans les nombreux tableaux qui suivent. 
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Dans les tables de corrélation, chacune des colonnes, portant en tête l’indi- 
calion d’un intervalle partiel de variation de d et, entre parenthèses, l'écart 
d— d, du centre de l'intervalle à la moyenne, fournit la répartition « liée » de 
l’autre variable a, b ou g. Au bas de la colonne est noté le total des observations 
de valeurs de d de l’intervalle partiel susdit et la moyenne liée de a, b ou g en d: 

Chacune des lignes donne les renseignements analogues, sauf à renverser le 
rôle des variables. 

Ces tableaux de calcul des écarts quadratiques ou dispersions se passent de 
commentaires, de même que les tableaux de calcul des écarts rectangulaires 
reproduisant, sous une forme différente, les tables de corrélation. 


REPRÉSENTATION GRAPHIQUE. — Les points expérimentaux, dont les coor- 
données sont définies par les moyennes liées et qui, théoriquement et de façon 
générale, sont situés sur les lignes de régression, ont été représentés géométri- 
quement sur des graphiques. On y a tracé les droites de régression et la droite 


idéale de la relation fonctionnelle linéaire qui correspondrait, dans chaque cas, 
àp=lt (Pl Tax) 


= 


né 7", cts 
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Gisement d’Angleur 
[Anthraconauta belgica Auctorum.] 


1. Tableau de corrélation entre les caractères d et &. 


a 
a (&— do) Fréquences | Moyennes 
G@— @o) | 43 | 85 | 5-7 | 7-9 | 9-41 | 11-43 | 13-15 | 15-17 47-19 a liées gen 
(4,20 | 2,0) | (-0,21) | (+1,89) | (43,89) | (+5,89) | (+7,89) | (:9,89) |(+11,89) 
1-3 
een AR TES 40 3,30 
3-5 
ie) 16e Se 4 1 86 5,07 
5-7 
dr 10 VAUT 1 67 7,91 
7-9 
Go 4 6 1 m 9,64 
9-11 
Gun 1 1 1 1 4 14,00 
11-43 : ta 
(1,37) 
13-15 
(@ 31) 1 1 16,00 
Requete us QUES ON COTES 2 1 2 1 209 Æ 
Moyennes © 
Ur de 54 4607 6,78 8 “ARE de 


We 
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\ ET 


2. Tableau de corrélation entre les caractères d et b. 


d 
| (d — do) Fréquences| Moyennes 
st | To | 92441 | 41-48 | 13-45 | 45-47 | 47-49 L Eu ii à 
(-0,41) | (+1,89) | (+3,89) | (+5,89) | (+7,89) | (+9,80) |(+14,89) 
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3. Tableau de corrélation entre les caractères & et g. 
d 
g (4— do) Fréquences | Moyennes 
G— 9) | 43 | 35 | 57 | vo | o14 | aus | 13-15 | 45-17 lav40 g liées d en g 
(-4,08) | (-2,08) | (-0.08) | (+1,92) | (13,92) | (45,92) | (27,92) | (49,92) |(+44,92) 
se 11 @WL A0 We 21 2,95 
(30) 44 4 5 1 51 4,12 
cL0180) 14 25 4 1 44 5,54 
(0:81) 8, tam 4 38 7,16 
( te ) 2 2 20 4 28 1,86 
ED 5. 6 11 9,09 
(3.60) TE UE 1 5 10,80 
(61) 1 1 12,00 
Gé) 1 1 16,00 
(6.61) RAP 2 17,00 
MUR Ts LD TON tes (He leur 2 1 2 1 202 
Moyennes gg 2151151718 5.068 15 04 7,50 10 10,50 


liées g en 4 


d3 
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1£‘6 68 L e 6 9 ë San 
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sauu9ÂO |SoouanPaa (0p —p) 7) 
p 


2. Tableau de corrélation entre les caractères d et 0. 
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Moyennes 
liées 4 en à 
,2 
4,10 
4,68 
6,00 
7,00 
7,87 
8,79 
9,71 
10,60 
11,43 
11,67 


13,00 


Fréquences 
b 
3 

10 
22 
20 
au 
30 
19 
7 
0 
# 
3 


149 


à 
(d — do) 
3 


11,67 


12 
1625 


10,36 


Ai 


14 
8,79 


8,14 


29 


24 
a, 21 


15 
5,87 


19 
5,32 


13 
4,54 


3,83 


12,5-13,5 
(5,44) 

413,5-14,5 
(6,44) 


(b 
2 
( 
3 
( 
4 
( 
5 
( 
6 
( 
7 
( 
8 


(1,44) 
d 


14,5-15,5 


Moyennes 
liées ben 4 


Fréquences 


BE) 
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$ 2. CALCUL DES MOYENNES,. 


Gisement de Wérister. 
[Anthraconauta minima HINp (non LUupwi1G), espèce de référence.] 


1 
dy A(d) = = PATP + 3 + 22 + 150 + 324 4 305 + 384 + 392 + 206 


295 -L 210 + 132 + 12 + 14] —6,31 


1 
Go = A (Q) = 25 (175 + 18 + 116 4 845 + 320 4 405 4- 300 147 +112 + 18 + 41] = 4,61 


1 
b= AG) = 5 (2175 + 2 + 20 + 120 + 340 + 400 + 330 + 420 + 344 


+ 225 + 170 + 44 + 156 + 14] — 6,44 


al 
go = A (9) = 2e LAS + 89 + 284 + 857 + 372 4 235 + 00 + 28 + 8] — 5,63 
499 
Gisement d’Angleur. 
[Anthraconauta belgica Auctorum.] 
bn . [28 + 288 -L 300 + 392 4 180 + 24414 + 32 P18] 6,11 
4 
de = A(@) = 35; [80 + 280 + 276 + 384 + 170 + 24 + 14 + 32 + 18] — 6,08 
G=A(G) = 2 80 + 344 + 402 + 88 + 40 4 14] — 4,63 
il 
bo À (D) = 35 (20 + 272 + 818 + 320 + 250 + 190 + 14 + 16 + 18] = 6,45 
1 


go = A (9) = 35 (101 + 21 + 102 + 182 + 152 + 140 + 66 + 35 + 8 + 9 + 20] — 3,89 


Gisement de Stud. 
[Anthraconauta belgica Auctorum.] 


dé = A (Q) = (18 + 52 + 05 + 00 + 168 + 282 + 126 + 140 + 182 4 86 + 80] — 7,3 

ao = A (a) = [T2 + 142 + 150 + 188 + 208 ++ 88 + 36] = 5,36 

b, = A (b) 1 [9 + 40 +410+ 120 119 + 240 + 171 + 70 + 410 + 84 + 39 + 15] = 7,56. 
ds À (0) = ra (26 + 90 + 120 + 2 25 + 114 + 63 + 24] — 4,44 


() Sur 209 individus. 
(2) et (*) Sur 202 individus. 


() 


(°) 
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8 3. GALCGUL DES DISPERSIONS. 


Gisement de Wérister. 
[Anthraconauta minima Hip (non LupwiG), espèce de référence. 


Calcul de 56, 


d | dæ | Fréquences F X d? 
1,5 2,25 3 6,75 
2,5 6,25 at 68,75 
3 9 412,25 50 612,50 
4,5 20,25 81 1.640,25 
5,5 30,25 73 2.208,25 
6,5 42,25 64 2.104,— 
7,5 56,25 56 3.150.— 
8,5 12,29 31 2.613,25 
955 90,25 25 2.256,25 
40/56 140,25 21 919,20 
MAPS 19225 12 1587, — 
1252810020 1 156,25 
1 54180720) — _ 
14,5 210,25 4 210,25 
435 D — 19.588,75 
A(d?) = 45,0316 
[A (d)]? = 39,8161 
Dre 100 
Er VER 
Calcul de 6, 
a | a? | Fréquences FX a? 
4,5 220 13 29,25 
a 6,25 58 362,50 
PA) 12:25 445 1.408,75 
4,5 20,25 80 1:620;— 
5,9 30,25 81 2.450,25 
6,5 42,25 50 24112750 
0,0) 56,25 21 1.181,25 
8,5 2525 14 41.011,50 
9,5 90,25 2 180,50 
10,5 110,25 — — 
LES 132,25 dl 132,25 
435 2 — 10.488,75 


A (a) — 241121 
[A (a)]? = 21,2521 
2 — 2,8600 

1,69 


a 
| 


a 
Il 


58 
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Calcul de 5, 


d | b? | Fréquences IR SQUE 
4,10 APS) 72 4,5 
2,5 6,25 10 62,5 
8,0 12725 40 490, -— 
4,5 20,25 85 1124525 
lojhis) 30,25 80 ArARUS= 
6,5 42,25 55 2 323,75 
Hal 60 31310200 
8,5 12525 43 3.106,75 
925 90,25 25 220010 
10/5004140729 117 1 874,25 
1155 0M62725 4 529; — 
12 50025 13 2.031,25 
43,5 182,25 — en 
14,5 : 210,25 4 210,25 
435 Z— 20.404,75 
A (b?) — 46,9126 
[A (b)]? = 41,4736 
où — D,4390 
Op — 2:99 
Calcul de 5, 
g | a | Fréquences F X.92 
s Ds) 2725 89 87,15 
2,0 6,25 EL 131,29 
AP) 19925 119 1.457,75 
4,5 20,25 93 1.883,25 
555 30,25 47 1.424,75 
6,5 42,25 15 633,75 
1,5 0056,25 4 225,00 
8,5 712725 1 72,25 
435 D — 6.512,75 . 
Ag?) = 14,9701 
CA (a)1=:13,1769 
CU 0e 


Il 


1,34 


DU HOUILLER DE LA BELGIQUE 


Gisement d’Angleur 
[Anthraconauta belgica Auctorum.] 


Calcul de 5, (1). 


U/ | d? | Fréquences Pod: 
2 4 14 56 
4 16 72 1452 
6 36 50 1.800 
8 64 49 3.136 
10 100 18 1.800 
12 144 2 288 

14 196 1 196 
16 256 2 512 

18 324 1 324 

209 X — 9.264 
ANOON= 218253 
PANIER 75321 
Gr 010082 
O4 — 2,64 
Calcul de 5, (2). 

d | æ Fréquences F X dæ 
2 4 15 60 
4 16 70 1.120 
6 36 46 1.656 
8 64 48 3.072 
10 100 17 1.700 
12 144 2 288 
14 196 ‘h 196 
16 256 2 512 
18 324 1 324 

202 2 — 8.928 
(d?) = 44,1980 
[A (d)]? = 36,9664 
EE SE 
da 00:00 


(1) Sur 209 individus. 
(2) Sur 202 individus. 


D9 
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Calcul de 5, 


a | a? | Fréquences FX a? 
2 + 40 160 
4 16 86 1.376 
6 36 67 2.412 
8 64 11 704 
10 100 4 400 
12 144 — — 
14 196 1 196 .. 

209 25.248 
A (a) — 25,1000 
[A (a)? ='21,4369 
CAS: 0681 
(A LC L 
Calcul de 6, 

b | b? | Fréquences F X 0? 
2 4 10 40 
4 16 68 1.088 
6 36 53 1.908 
8 64 40 2.560 
40 100 25 2.500 
12 144 10 1.440 
14 196 1 196 

16 256 al 256 

18 324 “l 324 

209 > — 4102912 
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Calcul de ‘y 


q | g | Fréquences F X g? 
4,5 2,25 21 47,25 
230 6,25 51 236,25 
08S 12,25 44 539 ,— 
4,5 20,25 38 769,50 
575 30,25 28 847 ,— 
6,5 42,25 44 464,75 
155 56,25 5 281,25 
8,5 12520 1 ASS) 
9,5 90,25 1 90,25 
10,5 110,25 2 220,50 

202 Z — 3.568,00 


A (g°) — 17,6633 
[A (g)]° = 15,1321 


a 
| 


2,5312 
sy = 1,59 


Gisement de Stud. 
[Anthraconauta belgica Auctorum.] 


Calcul de 6, 


di | d? | Fréquences FX? 
3 9 6 54 
4 16 13 208 
(5) 25 19 475 
6 36 15 540 
ÿ 49 24 1.176 
8 64 29 1.856 
9 81 14 41.134 

10 100 14 1,100 

44 121 12 1552 
12 144 3 432 

13 169 3 507 

149 2 — 8.934 
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Calcul de 5, 


a | a? | Fréquences F x d? 
3 9 24 216 
4 16 28 448 
Hi) 25 30 750 
6 36 23 828 
7 49 29 1.421 
8 64 11 704 
9 81 4 324 | 

149 = 4.691 
AT) = 19187882 
[A (a)]? = 28,7296 
C2 1086 
60 O0 
Calcul de 5, 

b | b? | Fréquences F x 0? 
3 9 3 ot 
4 16 10 160 
o 25 22 550 
6 36 20 720 
Fi 49 47 833 
8 64 30 1.920 


© 
[ee] 
= 
2 
© 
= 

ot 
[0] 

de] 


= 
D = 
= à 
æ 4 
= + 
pe 
1 © 
= > 
8£ 
œ@ © 


13 169 6) 507 

14 196 ee = 

15 225 { 225 
149 Z— 9.399 


UT = 
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Calcul de Sy 


q | g? | Fréquences F X g° 
2 4 13 52 
3 9 30 270 
4 16 30 480 
5 25 45 4,425 
6 36 19 684 
t - 49 9 441 
8 64 3 192 

149 È — 3.244 


A (g°) = 21,7718 


[A (9)]? = 19,7136 
Gb NUE 
Do LAS 


8 4. — CALCUL DES ÉCARTS RECTANGULAIRES ET DES COEFFICIENTS 
DE CORRÉLATION ET DE RÉGRESSION. 


Gisement de Wérister. 
[Anthraconauta minima Hip (non LUDWwIG), espèce de référence.] 


Calcul de l’écart rectangulaire E?7, et des coefficients p4, et de régression 
de «a en d et d en «. 


d—& a— &o | (4 — do) (a — &y) | Fréquences F X (4 — do) (a — &o) 
—4,81 ll 14,9591 3 44,8773 
3 si 11,8491 9 106,6419 
3,81 oil 8,0391 2 16,0782 
Xl ol 8,7391 1 8,7391 
ol ten 5,9291 40 237,1640 
4 sai] 3,1191 9 28,0719 
il oil 3,8191 15 57.2865 
—1,81 Et 2,0091 60 120,5460 
el OU 0,1991 6 1,1946 
—0,81 fl 1,7091 1 1,7091 
—0,81 NE 0,8991 40 35,9640 
Der = 0,0891 30 2,6730 
—0,81 +0,89 —0,7209 2 —1,4418 
+0,19 ti —0,2109 5 —1,0545 
+0,19 in) —0,0209 33 —0,6897 
+0,19 0,89 0,1691 25 4,2275 

ATODOTIEL.-. 22e 281 


2 
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Reporter: 281 
+0,19 1,89 0,3592 1 0,3591 
1,19 nl —1,3209 1 —1,3209 
1,19 bi —0,1309 11 —1,4399 
1,19 0,89 1,0591 34 36.0094 
1,19 1,89 2,2491 10 22,1910 
2,19 0,89 1,9491 18 35,0838 
2,19 1,89 41391 19 78.6429 
3,19 0,89 2,8391 o DEN 5 6782 
3,19 1,89 6,0291 Fame 72.3492 
3,19 2,89 9,2191 10 92,1910 
3,19 3,89 12,4091 ñ 12,4091 
419 1,89 7,9191 ü 55,4337 
4,19 2,89 121091 10 121,0910 
4,19 3,89 16.2991 4 65,1964 
5,19 1,89 9,8091 1 9,8091 
5,19 2,89 14,9991 1 14,9991 
5,19 3,89 20,1891 9 181,7019 
5,19 4,89 25,3791 1 25,3791 
6,19 4,89 30,2691 1 30,2691 
8,19 6,89 56,4291 1 56,4291 
435 S 2157472485 
E2, — SE » — 3,6201 
PT 


ue a. Eôs 3,6201 
Coefficient de corrélation : pau — = —— —= 0,940 
Sa Su 3,8D92 
Coefficient de régression d papes 0 2 6694 
e : régression de & en d: — =: — 0,69 
oefficient de rég 0 Paa 5,  B218S : 
Coefficient de régression de d QE 
oefficient de régression de d en &: — 1,262 = 0,7 
: Ë 58600 
NS : c, 169 
Coëfficient angulaire pour = 1 : == — — 0,741 
gau 228 
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Calcul de l’écart rectangulaire E2y et des coefficients par et de régression 
de b en d et d en b. 


d — do b— Do | (4 — do) (db — bo) | Fréquences FX (4 — do) (b — do) 
hat HS 04 23,7614 2 7,5228 
ait ITE5:97 18,9514 1 18,9514 
381 inc 8704 15,0114 7 105,0798 
381 O0: 11,2014 4 44,8056 
2,81, DIE 08 11,0714 1 11,0714 
—2,81 d11029,94 8,2614 33 272,6262 
sit 19% 5,4514 16 87,2224 
= Re 71314 1 7,1314 
LE 0) 5,3214 3 15,9642 
ER A 3,514 56 196,6384 
48h 09% 1,7014 e1 35,7294 
Ha 19% 1,5714 12 18,8568 
EDS dx 0,7614 45 34,2630 
0,81 +0,06 —0,0486 16 —0,7776 
+019 | 19% —0,3686 1 —0,3686 
+0,19  —0,94 —0,1786 11 —1,9646 
+019 +0,06 0,014 32 0,3648 
+0,19 +1,06 0,2014 18 3,6252 
+019 2,06 0,3914 2 0,7828 
41,19  —0,9% —1,1186 2 22970 
+1,19 +0,06 0,0714 6 0,4284 
1.19 1,06 1,2614 32 10,3648 
1,19 2,06 2,4514 16 39,2224 
2,19 - — 0.94 2 —2,0586 1 —2,0586 
219 +0,06 0,1314 1 0,1314 
2,19 1,06 2,3214 10 23,2140 
2,19 2,06 4,5114 18 81,2052 
2,19 3,06 6,7014 6 40,2084 
2,19 4,06 8,8914 1 8,8914 
3,19 2,06 6,5714 7 15,9998 
3,19 3,06 9,7614 12 1171368 
3,19 4,06 12,9514 6 77:7084 


A reporter. … … 400 
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FRÉDOL EME 400 
4,19 3,06 12,8214 Y 89,7498 
4,19 4,06 17,0114 9 153,1026 
4,19 5,06 21,2014 3 63,6042 
4,19 6,06 25,3914 2 50,7828 
5,19 4,06 21,0714 1 21,0714 
5,19 5,06 26,2614 1 26,2614 
5,19 6,06 31,4514 10 314,510 
6,19 6,06 37,5114 1 37,5114 
8,19 8,06 66,0114 1 66,0114 
435 Y =2190,3490 
E, — Re 50809 
sh 466 
y Op 0 2,99 —19,9184 
Coefficient de corrélation : puy = — 0,948 
Sa So 
, , CA 5,0353 
Coefticient de régression de à en d: — — 0,965 
| P& 52158 
. "0 : Cu 5,0353 
Coefficient de régression de d en b:  — — 0,925 — 
: Fe ne 5.430000 
£ L 67 L283 
Coefficient angulaire pour p = 1 : — — 1,022 
Sa 2,28 
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Calcul de l’écart rectangulaire E?3, et des coefficients p4, et de régression 


de g en d et d en g. 


a 
| 
d — do g — Go | (4 — do) (g — go) | Fréquences EX (4 — do) (g — go) 
4,81 2913 10,2453 3 30,7359 
—3,81 213 8,1153 11 89,2683 
ie 213 5,9853 24 143,6472 
PH = 3,1753 25 79,3825 
—?,81 —0,13 0,3653 1 0,3653 
el EL 3 3,8553 3,8553 
—1,81 18 2,0453 64 130,8992 
—1,81 —0,13 0,2353 10 2,3530 
—1,81 +0,87 —1,5747 6 —9,4482 
—0,81 Ne 0,9153 27 24,7131 
0/81 2UNT-20;13 0,1053 & 4,1226 
—0,81 +0,87 —0,7047 4 —2,8188 
+0,19 ‘1,13 —0,2147 20,214? 
0,19 : 0,13 —0,0247 45 11115 
0,19 +0,87 0,1653 18 2,9754 
1519 2000-20 18 —0,1547 19 —2,9393 
1,19 0,87 1,0353 33 34,1649 
1,19 1,87 2,8253 8,9012 
219 SENS ;13 —0,2847 —0,5647 
2,19 0,87 1,9053 26 19,5378 
2.19 1,87 4,0953 36,8577 
3,19 0,87 2,7753 11,1012 
3,19 1,87 5,9653 17 101,4101 
3,19 2,87 9,1553 4 36,6212 
4,19 0,87 3,6453 8 7,2906 
1,19 1,87 7,8353 13 101,8589 
4,19 2,87 12,0253 6 72,1518 
5,19 1,87 9,7053 ! 38,8212 
5,19 2,87 14,8953 5 74,4765 
5,19 3,87 20,0853 3 60,2559 
6,19 3,87 23,9553 1 23,9553 
8,19 4,87 39,8853 1 39,8853 
435 Z=1192,8102 
Eôg — 18 — 2,7421 
Sa Sy = 2.28 X 1,34 — 3,0552 
Coeflicient d rélation : buy — Leg nuot Midi 0,898 
oelicien e CorrelaliOn : Pag 75, 3,0552 be 
k . Rue ei — 
Coefficient de régression de g en d'! Pa = — pose d 26 
. . ra RTE E 
Coefficient de régression de d en g : 84 % — 17932 — 1,529 = 0,654 
Coefficient angulaire pour p = 1: RSR 0,588 
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Gisement d’Angleur. 
[Anthraconauta belgica Auctorum.| 


Calcul de l’écart rectangulaire E?4, et des coefficients pq, et de régression 
de «a en d et à en a. 


4 — do a — dy (4 — do) (a — &o) | Fréquences FX (4— do) (a — &o) 
—4,11 —?,63 10,8095 14 --151,3302 
—2,1l —?,63 5,5493 —26 5 144,2818 
ill —0,63 1,3293 46 61,1478 
in —0,63 0,0693 35 2,4255 
—0,11 +1,37 —0,1507 15 —2,2605 

1,89 —0,63 —1,1907 4 —4,7628 

1,89 31,37 2,5893 41 106,1613 

1,89 So 6,3693 3 4 25,4778 

3,89 —0,63 —?,4507 1 —2,4507 

3,89 1,37 5,3293 10 53,2930 

3,89 3,00 13,1093 6 78,6558 

3,89 Bron 20,8893 1 20,8893 

5,89 1,37 8,0693 1 8,0693 

5,89 5,37 31,6293 1 31,6293 

7,89 ad 26,5893 1 26,5893 

9,89 Hat 53,1093 1 53,1093 

9,89 9,37 92,6693 1 92,6693 

11,89 eye 63,8493 d 63,8493 | 
209 Y=910,1037 | 
Da = 4,3946 
Sy Sa = 2,64 X 1,91 —5,0424 
at Lu 200046 

Coefticient de corrélation : pgx — — 0,864 


Sy Sy -0,0424 


3546 
Coeflicient de régression de & en d: re — —10:02S 
; $ D £ = A — A) 

L fé pire 


Coeficient de régression de d LR A 
oelicien e regresslon de enr: — = = 1,109 = —— 
# FPS 36631 0,841 


. . ÿ, 
Coefficient angulaire pour p = 1 : = — 0,723 
Sa 
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Calcul de l’écart rectangulaire E27, et des coefficients par et de régression 
de b en d et d en b. 


d — do D— (d — do) (bd — bo) | Fréquences EF X (4 — do) (b — Do) 
QU 18,2895 7 128,0265 | 
PME 7 10,0695 7 70,4865 
DS PU 35 9,3895 3 28,1685 
Ep tie ES 5,1695 59 305,0005 
ob 35 0,9495 10 9,4950 
FM 0,2695 2 0,5390 | 
= INT ds 0,0495 38 1,8810 
0,11 1,55 __0,1705 10 _1,7050 
1 80 POS n 45 —_0,8505 5 Dy.9525 
1,89 1,55 2,9295 28 82,0260 
1,89 ne 6,7095 16 107,3520 
3,89 Ds 6,0295 2 12,0590 
3,89 3,55 13,8095 9 124,2855 
3,89 5,55 21 5895 7 1511265 
5,89 8,55 32,6895 2 65,3790 
7,89 $ ss 13,7895 1 13,7895 
9,89 7,55 72,6695 1 72,6695 
9,89 Â1:55 1142295 1 1142295 
11.89 9,55 113,5495 1 113,5495 
209 S 1426,1055 
Er —= 6,8235 
Su = 2,64 X 2,78 = 7,3392 
o PP ON cost 
Coefïicient de corrélation : bg — 7.3302 —°0,93 
cond end», = PE Gore 
OefliCieni de récression de Où CNP AE ble mamie Un 6) 
RE RENNES 
É hient deré Lo An ELU, CSL PP SN © 1 
oelicien erré2TressiIOn Aie arenr0 : — = M OC PA 
8 Pdè D TT372 1,134 
Coefficient angulai MR RER ee 
Oeilicient angulaire pour 0 — 5 = — = Do 
8 RO emo 
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Calcul de l’écart rectangulaire E, et des coefficients p4 et de régression 


de g en d et d en g. 


d— do | g—go | (4 — do)(g — go) | Fréquences F X (4— do) (g — go) 
2,08 2839 9,7512 11 107,2632 
—4,08 —1,39 5,6712 4 22,6848 
—,08 —?,39 4,9712 10 49,7120 
—2,08 —1,39 2,8912 41 118,539 
—2,08 —0,39 0,8112 14 — 11,3568 
—2,08 0,61 —1,2688 Rare —3,8064 
—?,08 +1,61 —3,3188 8 —6,6976 
—0,08 —1,39 0,1112 5 0,5560 
—-0,08 —0,39 0,0312 85 0,7800 
—0,08 0,61 —0,5488 14 —(0,6832 
—0,08 1,61 —0,1288 ê —0,2576 
1,92 —1,39 —?,6688 1 —?,6688 
1:92 —(0,39 —0,7488 4 —?,9952 
1:92 0,61 1,1712 17 19,9104 
1,92 1,61 3,0912 20 61,8240 
1,92 2,61 5,0112 5 25,0560 
LE 3,61 6,9312 1 6,9312 
3,92 —(0,39 —1,5288 1 —1,5288 
3,92 0,61 2,3912 4 9,5648 
3,92 1,61 6,3112 4 25,2448 
3:02 2,61 10,2312 6 61,3872 
3,92 3,61 14,1512 8 28,3024 
5,92 3,61 e1,3712 1 21,3712 
5,92 4,61 27,8918 1 278912 
7298 3,61 28,912 1 28,5912 
9:08 5,61 55,6512 1 55,6512 
9,92 6,61 65,5712 il 65,5712 
11,92 6,61 78,7912 3} 78,7912 
202 X = 807,7424 
as = 3.9987 
Sa Sy — 2,69 X 1,59 — 4,2771 
. FA 3,9987 
Coefficient de corrélation : 949 — ITA — 0,935 
La] 
Coefticient de régression de g en d : Pay e = =. — 0,553 
IAE .. ‘og 80987 1 
Coefticient de régression de 4 en g : pag — = SE — 1,580 — 0,633 


Sg 


Coefficient angulaire pour p = 1: en 0,591 
a 
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Gisement de Stud. 
[Anthraconauta belgica Auctorum.] 


Calcul de l’écart rectangulaire L?,, et des coefficients pau et de régression 
de «a en d et d en a. 


A — dy a— ao | (d— do) (a — &o) | Fréquences | EF X (4 — do) (a — &o) 
—4,37 06 10,3132 6 61,8792 
337 26 7,9532 11 87,4852 
—3,37 —1,36 4,5832 2 9,1664 
—2,37 —? ,36 5,5982 5 27,9660 
Cl —1,36 3,2232 12 38,6784 
27 —0,36 0,8532 2 1,7064 
—1,37 —?,36 3,2332 2 6,4664 
137 —1,36 1,8632 9 16,7688 
tot —0,36 0,4932 3 1,4796 
hr +0,64 —0,8768 1 —0,8768 
—-0,37 —1,36 0,5032 n 2,0128 
—0,37 —0,36 0,1332 10 1,3320 
—0,37 +0,64 —0,2368 8 —1,8944 
0:37 1,64 —0,6068 Û —1,2136 
+0,63 —1,36 —0,8568 L —0,8568 
+0,63 —0,36 —0,2268 14 —3,1752 
0,63 0,64 0,4032 8 3,2256 
0,63 1,64 1,0332 6 6,1992 
1,63 —0,36 —0,5868 Al —0,5868 
1,63 0,64 1,0432 4 4,1728 
1,63 1,64 2,673 9 24,0588 
2,63 0,64 1,6832 2 3,3664 
2,63 1,64 4,313? 5 21,5660 
2,63 2,64 6,9432 ui 27,11788 
3,63 1,64 5,9532 ? 11,6724 
3,63 2,64 9,5832 3 28,7496 
3,63 3,64 13,8132 8 26,4264 
4,63 __ €,64 12,2232 è 24,4464 
4,63 3,64 16,8532 1 16,8532 
5,63 2,64 14,8632 8 29,7264 
5,63 3,64 20,4932 1 10,9864 
149 2 = 545,600 
E2, —9,6615 
00, 2,010C1,00119,0842 
euh 2-7 00088 à ne 3.6614 bé 
Coefficient de corrélation : pau = _n oi 3 
« : MS rs ; 
Coefficient de 70. de «& en d: Pau ss = 56428 — 0,649 
, ee Tee 0010 îl 
Coefficient de régression de d en & : Pan — — Et 1,330 — = 
CR 100 0,792 
Goefficient angulaire pour bp —1: FEMMES 0,700 


Sa 2,37 
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Calcul de l'écart rectangulaire E°4, et des coefficients pa et de régression 
de b en dd et d en b. 


d — dy | 9 — go | (&— do)(g — go) | Fréquences F X (4 — do) (g — go) 
ee 19,9272 2 39,8544 
LAS CO OËSE 15,5572 3 16,6716 
STE 66 11,1872 1 11,1872 
2337 gte56 15,3672 1 77745,3672 
"5 97 MESSE, 11,9972 TRE | 47,9888 
23,370 MMED:SS 8,6272 8 69,0176 
29878" DE 8,4372 2 16,8744 
5870, NE 06 6,0672 10 60,6720 
297 aa 3,6972 6 22,1832 
2,37, 26180 56 1.3272 1 1,3272 
an MEL 56 4,877 1 1,877 
Sr Me 3,5072 3 10,5216 
LATE MEAUSE 2,1372 9 19,2348 
4:87 wer2056 __ 0,62 1 0,7672 
21,570 ND 22 —0,6028 1 —0,6028 
0,37 NOR 56 0,572 4 2,3088 
me") 0,2072 12 2,4864 
087 0 —0,1628 7 —1,1396 
ONLINE —0,5328 1 —0,5328 
+0,63 21,56 —0,9828 1 —0,9828 
+0,63 —0,56 —0,3528 3 —1,0584 
+0,63 0,44 0,2772 17 4,7124 
0,64 1,44 0,9072 7 6,3504 
0,63 2,44 1,5372 À. 1,5372 
1,63 0,44 0,7172 5 3,5860 
1,63 1,44 2,3472 ? 16,4304 | 
1,63 2,44 3,9772 2 7,9544 | 
2,63 1,44 3,7872 3 11,3616 | 
2,63 2,44 6,4172 2 12,8344 | | 
2,63 3,44 9,0472 5 45,2360 
2,63 4,4% 11,6772 1 11,6772 | 
3,63 1,44 5,2272 1 5,2272 
3,63 2,44 8,8572 2 17,144 
ESS ll 
ANTODONEL, ra 134 | 
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9 ine/0oitisco 280 19 | grieium 


HÉDOR- A. E. 134 
RAS 
16,171 
19,7472 
_ 15,927 
20,552 
22,9972 
30,6272 
11,8872 


19,988 
48,3516 
39,944 
15,9272 
a1,1144 
22,9972 EX 
30,6272 
11,8872 


DE 0 col ro cod Ke 


149 © — 803,902 


ir #4 EE é ’ 
a dés fu: Æ à == 5,3050 
1 HR cu 237 X 2,43 53 59 


0 rétati s _ Ban 53050 0,937 paie 
orré ion : Par a auO + 


, ÿ' C1) - 


ce 5,3959 


it. er ea. 1 


Sa 
P 
1 
#g = : 
: ( = 
| Hi Î 
1 
: nc 
’ Î L 
(1 l 
LI 
| 
0: 
t,L1n7 
1 
x 
me , { M UE ) 41 
ue : Vy ta *» aù FAO 1 Lite ) 
ve , L 
v L } pe ; } 1e FA 
0 CO ED 08 NO #YTNS Ab Jus EN 
= : 


= +} = 0 10004 s1isirgus Ipsioiso 
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Calcul de l’écart rectangulaire E?7, et des coefficients pag et de régression 
de 4 en d et d en g. 


dd db bo | dd) 0) | Fréquences | EF X (4 — do) (a — &o) 
| 
ST AE 6,2928 1 6,2928 
=87 ‘eu 10,6628 5 54,3140 
2957 MP LL 8,2228 5 41,140 
3,370, dx 4,8528 ? 33,9696 
3.370000 042 1,4828 1 1,488 
OST ULE 5,7828 2 #1,5656 
Hein ARE 3,4128 13 44,3664 
D NT RENE 1,0428 4 21712 
ASC EUr 3,3428 1 3,3428 
LH 207 RUE 1,9728 7 13,8096 
HT ONE CE 0,6028 6 3,6168 
Zi "+056 —0,7672 1 —0,7672 
2637 NE 0,5328 1 0,5328 
2087 "SN 0,1628 13 2,1164 
0,37 +0,56 —0,2072 8 —1,6576 
et —0,5772 2 1,154 
10,68 > —0,9072 1 —0,9072 
+0,63  —0,44 —0,2772 6 —1,6632 
0,63 0,56 0,3528 20 7,0560 
0,63 1,56 0,9828 2 1,9656 
1,63 0,56 0,9128 11 10,0408 
1,64 1,56 2,5428 3 7,6284 
2,63 0,56 1,4728 3 4,184 
2,63 1,56 4,1028 5 20,5140 
2,63 2,56 6,7328 3 20,1984 
3,63 0,56 2,038 2 4,0656 
3,63 1,56 5,6628 5 28,3140 
3,63 2,56 9,2928 4 37,112 
3,63 3,56 12,9228 1 12,9228 
4,63 1,56 7,2228 2 14,4456 
1,63 3,56 16,4828 FAT 16,4828 
5,63 2,56 14,4128 2 28,8256 
5,63 3,56 20,0428 d 20,0428 
149 5 = 447,632 
E, — 3,0043 
Coefficient de corrélation : oy = äo = — 0,882 
À é PAT eee DSTI di 
Coefticient de régression d PU ere 
10601CI1en e regression «qe { en : Pag ea = 56428 —= V,99£ 
sa Eu Sa  3,0043 | 1 
2oefficient de régression de d en g : Pay e = 5058 — 1,460 — 0,685 
Q 
Coefticient angulaire pour 9 = 1 : U — = 0,603 
Ta 201 


oi 
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CONCLUSIONS. 
\ 


Il a été réuni, dans le tableau ci-dessous, les valeurs des coefficients de 
corrélation pu; Om Lay Pour chacune des trois populations étudiées. 


ne Coeflicients | 
Pda db Pag 
A.— Wérister 0,940 0,948 0,898 
A.— Angleur 0,864 0,930 0,935 
A Stud 0,931 0,937 0,886 


Ensuite il a été calculé les écarts quadratiques dont l’utilisation serait légi- 
time si les distributions étaient régies par une loi de Bravars. Rappelons que ces 
écarts quadratiques, calculés au moyen de la formule 


nous fixent l’ordre de grandeur des écarts à partir de la valeur théorique exacte 


de p. 


Ecarts quadratiques 
Origine 
Spa S pa S Pag 
A.— Wérister 0,0058 0,0049 0,0093 
A.— Angleur 0,0175 0.0093 0,0087 
A.— Stud 0,0109 0,0099 0,0176 


Remarque. — Dans la présente étude interviennent simultanément plusieurs 
caractères : d, a, b, g. PEaRsoN (*) a fait observer qu'une telle étude nécessite, 
en général, l'introduction d'autres coefficients que les coefficients de corrélation 
relatifs aux combinaisons des divers caractères, pris deux à deux [G. U. Yuze 
et M. Kenpazz (°*)]. 

. Maïs la quasi-linéarité des relations entre les couples de variables correspon- 
dantes suffit ici pour établir les conclusions relatives à la stabilité de forme 
géométrique des coquilles, qu’on trouvera plus loin. 


(3) Unny YuLE, G., et KenDaL, M. G., 1945, pp. 261-286. 
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E. — ETUDE DES INDICES CARACTERISTIQUES. 
(Rapport entre mesures des caractères des individus. Distribution de fréquences.) 


Méruone. — Selon les méthodes d'établissement de relevés statistiques por- 


tant sur des variables continues, exposés au chapitre IL, il a été procédé à un 
relevé direct des distributions de fréquences des indices caractéristiques. 


Ces indices caractéristiques, étudiés dans les trois cas pour À. — Angleur, 
A. — Stud et A. — Wérister, sont les rapports 
7 d ia b 
IUO AOÛ, 100 X —, 
: a 9 


obtenus en effectuant les quotients de deux des grandeurs caractéristiques 
mesurées sur les coquilles. 

Désignant par æ, de façon générale, ces indices, il a été compté dans chaque 
intervalle partiel de longueur h convenablement choisie, centré sur les points x, 
les fréquences absolues o(x;)h; ce qui a permis de déterminer les valeurs expé- 
rimentales de la fonction à paliers o(x,). 

Les moyennes et dispersions permettent de déterminer la forme analytique 
de la loi de Gauss : 


qui, en première et suffisante approximation, réalise l'ajustement de la distri- 
bution des fréquences absolues. 

Pour chacune des trois populations étudiées et pour chacun des indices 
caractéristiques, les résultats (relevés bruts), les calculs et ajustements sont 
présentés ci-dessous de façon uniforme. 

N est le nombre de coquilles mesurées, n le nombre d’intervalles partiels et 
R la longueur de ces intervalles, À et ç« la moyenne et la dispersion. Le coeffi- 


. . ç + 17 + . . . 
cient de concentration + donne une idée de l'importance relative de la dispersion 


par rapport à la moyenne. 
. 


est l’ordonnée maximum de la courbe de Gauss qui réalise l’ajuste- 


SV2r 
ment des distributions représentées dans les figures des planches X à XIV sous 
la forme de polygones de rectangles, Ce maximum théorique est atteint pour 
une valeur de æ égale à A. 

Dans les tableaux rassemblant les relevés et les résultats des ajustements, 
les titres des colonnes ont la signification suivante : 


1) æ,: centres des intervalles partiels de longueur h; 


2) X;: valeurs correspondantes de la variable auxiliaire X définie par 
x=hX EE; 


DU HOUILLER DE LA BELGIQUE 71 


3) o(x,)h : fréquences absolues; 
4) o(x,) : ordonnées de la fonction à paliers; 
5) æ,— A : écarts de la moyenne; 
%; — À : den eZ ; ; 
6) - seconde variable auxiliaire utilisée pour l'emploi des tables 


donnant les valeurs de la fonction y — Wie, # 
T 


lè 


7) ——e 2 : valeurs tirées d’une table; 
V2 
lè 
8) - —e : 2 : valeurs ajustées de la fonction + (x) ordonnées de la 
5 V2r 


courbe de Gauss, correspondant aux abscisses +. 


$ 5. DISTRIBUTION DE FRÉQUENCES. 


Gisement de Wérister. 
[Anthraconauta minima Hip (non LupWwiG), espèce de référence.] 


1. Indice 100 x 1 
N 
N — 433 Ae=t71,04 "11243 
(ox 
RMC oc — 10,76 
ç N 
h=5 Le 0,150 — — 16,054 
o 2) 
; € 42 
100 |. Xi o (x) A | o (@i) A UP home 
b | 1 ex 24 
50 — 5 1 —21,94 —2,04 0,04999 2,012 
55 —3 15 3 —16,94 —1,57 0,11650 4,688 
60 =? 51 10,2 —11,94 —1,11 0,21548 8,671 
65 —1 94 18,8 — 6,94 —0,64 0,32483 13,072 
70 0 96 19,2 — 1,94 —(,18 0,39222 15,784 
75 1 64% 12,8 3,06 0,28 0,38332 15,426 
80 2 36 re 8,06 0,75 0,30097 LA TULe 
35 3 34 6,8 13,06 1,21 0,19191 2180 
90 n 19 3,8 18,06 1,68 0,09742 3,920 
95 5 8 1,6 23,06 2,14 0,04058 1,633 
100 6 8 1,6 28,06 2,61 0,01326 0,533 
«| 105 7 il 0,2 33,06 3,07 0,00362 0,146 
110 8 = -— — — — — 
115 9 il 0.2 43,06 4,00 0,00013 0,005 
120 10 — EE = EE RE + 
125 11 1e 0,8 53,06 4,93 — —— 


Écart probable : 0,6745 x10,76= 7,258. 
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N=0453 A AS : 22102022 

ME 20 os = 13,58 

he T 0,099 patnbeg 

A cs V 27 
2 9 
100 = ARS o (æi) À p (ti) %i — A FAOTÉS TEE _ N._ (A 
4 TT 6 É V'2r V2 
90 10 1 0,2 gs 97 23,2 0,00116 0,038 
95 _9 0 0 40,77 3,05 0,00385 0,125 

100 8 1 0.2 245,7 2,67 0,01137 0,368 
105 de £ 0.8 —30,77 —_2,30 0,02833 0,917 
110 6 6 1,2 Les 7 1,93 0,06213 2,011 
115 Le 22 44 20,77 21,55 0,12022 3,891 
120 vf 34 6,8 2045, 97 4,18 0,19892 6,137 
125 40 8,0 __10,77 —_0,80 0,28969 9,375 
130 Lie 50 10,0 25,77 0,33 0,36341 11,761 
135 Mine “0 25 10,6 Lot —_0,06 0,39775 12,872 
140 0 78 15,6 1,23 0,32 0,37877 12,258 
145 1 47 9,4 9,23 0,69 0,31015 10,037 
150 2 58 6 14,23 1,06 0,22750 7,362 
155 3 16 3,2 19,23 1,44 0,14165 2,584 
160 n 7 1,4 24,23 1,81 0,07762 2,512 
165 5 6 1,2 29,23 2,18 0,03717 1,203 
170 6 6 1.2 34,23 2,56 0,01516 0,491 
175 7 3 0,6 39,23 2,93 0,00549 0,178 
180 8 0 0 14,23 3,31 0,00167 0,054 
185 9 1 0,2 19,23 3,68 0,00046 0,015 


Écart probable : 0,6745 x13,38=— 9,085. 
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; a 
3. Indice 100 x - 
% N 
NI= 2425 A=—1179:52 20722 
© 
n = 14 o = 20,51 
Ci N 
h = 10 F 10114 - VE: — 8,267 
lt? 
Li — À AMEL AT N le 
100 — = x; X: o (ui) ? (@i) Er A = nl — € 
" V’2r V2» 
115 6 7 0,7 —64,52 3,15 0,00283 0,059 
125 —5 3 0,3 —54,52 —2,66 0,01168 0,042 
135 =; 4 0,4 —44,52 1? 0,03802 0,788 
145 Æ 15 1,5 —34,52 68 0,09742 1,999 
155 _ 25 2,5 24,52 —1,20 0,19419 4,024 
165 Bar 68 6,8 —14,52 ar 0,31000 6,424 
175 0 93 9,3 165 0,2? 0,38911 8,063 
185 1 91 9,1 + 5,48 0,27 0,38429 7,963 
195 2 62 6,2 15,48 0,75 0,30097 6,237 
205 3 29 2,9 25,48 1,24 0,18506 3,835 
215 4 17 L7 35,48 1,73 0,08952 1,855 
295 5 9 0,9 45,48 2,29 0,03404 0,705 
235 6 1 0,1 55,48 27h 0,01017 0,211 
245 7 1 0,1 65,48 3,19 0,00229 0,047 


Écart probable : 0,6745 x 20,51 =13,833. 
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Gisement d’Angleur. 
[Anthraconauta belgica Auctorum.] 
É a 
1. Indice 100 x =. 
b 
N 
N — 209 À — 72,92 — — 14,8227 
(y 
h = 10 o — 14,10 
ç N 
= A0 —— 0,1934 = — 5,913 
A s V2 T 
li t 
« Gi — À 1 Re à N re 
1007 —a | Xi o(ai)h w(æi) otage eg -| 60 
| da Vox sV/ 2x 
50 —? 15 15 —22,92 1,6 0,11092 1,64 
60 —1 51 5,1 —12,92 —0,9 0,26609 3,94 
7 0 64 6,4 — 9,92 —0,2 0,39104 5,80 
80 1 40 4,0 7,08 0,5 0,35207 5:08 
90 2 27 Pa 7,08 ‘1 0,19419 2,88 
100 3 fa 6,7 27,08 1,9 0,06562 0,97 
110 4 2 0,2 37,08 2,6 0,01358 0,20 
120 5 0 0 47,08 39 0,00172 0,0? 
130 6 2 0,2 57,08 4 0,00013 — 
140 7 1 0,1 67,08 4,7 0,00001 — 
Écart probable : 0.6745 x14,10—9,51. 
; d 
2.4Indice 100% - 
a 
N — 209 À — 131,26 N = 9,9335 
h = 10 CL 0 
C2 N 
n —= 1413 — = (),1603 —— — 3,9629 
À C) V2 T 
2 dl? 
Xi — À 4 CRE N ET 
100-—axi|  X; o(æ:)h w (æ5) 2 — Al ——e SORA 
a 2 Von CJ V/2r 
77 —5 1 0,1 —54,26 —?2,6 0,01358 0,13 
87 —4 d 0,1 —44,26 —?,1 0,04398 0,44 
97 —3 5 0,5 —34,26 —1,6 0,11092 1,10 
107 =? 23 2,3 —24,26 11 0,20602 2,05 
117 —1 40 4,0 —14,26 —(0,7 0,31225 3,10 
127 0 48 4,8 4,26 —0,2 0,39104 3,88 
137 1 44 4,4 6,74 0,3 0,38139 3,79 
147 8 20 2,0 16,74 0,8 0,28969 2,88 
157 3 12 17e 26,74 13 0,17137 1,70 
167 4 5 0,5 36,74 17 0,09405 0,93 
177 5 3 0,3 46,74 22 0,03547 0,35 
187 6 ti 05 56,74 Qt 0,01042 0,10 
197 7 2 0,2 66,74 D 0,00238 0,02 


Écart probable : 0,6745 x21,04=— 14,19. 


DU HOUILLER DE LA BELGIQUE 81 


RU 100». 


(2 
NE 01 A 109,21 : — 7,1308 
—10 gs == 26,0 
n = 18 = 01573 AL el 2 
os V2r 
£ 2 
100 en Qi X: wo (%)R © (Œi) œ ee En _ EN. : rie EE . 

g - | k “ T (J Ve T 
80 =8 1 0,1 85,27 + 0,00172 0,01 

90 — 2 0,2 75,27 —2,9 0,00595 0,05 
190 — 1 0,1 —65,27 2,5 0,01753 0,14 
110 _& 3 0,3 —55,27 2,1 0,04398 0,34 
120 y 5 0,5 —45,27 Li 0,09405 0,73 
130 2 9 0,9 35,27 1,35 0,16055 1,24 
140 _ 13 1,3 25,27 1,0 0,24191 | 1,87 
150 Si 27 2,7 15,27 —0,6 0,33322 2,58 
160 0 12 42 — 5,27 —0,2 0,39104 3,02 
170 1 37 a 4,73 +02 0,39104 3,02 
180 2 18 1.8 14,73 +0,55 0,34264 2,65 
190 3 19 1,9 24,73 0,95 0,25403 1,96 
200 4 11 di 34,73 1,3 0,17137 1,32 
210 5 8 0,8 44,73 1,7 0,09405 0,73 
220 6 2 0,2 54,73 2,1 0,04398 0,34 
230 7 2 0,2 64,73 2,5 0,01753 0,14 
240 8 E ni 74,73 2,9 0,00595 0,05 
250 9 1 sn l 84,73 33 0,00172 0,01 


Écart probable : 0,6745 » 26—17.54. 
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Gisement de Stud. 
[Anthraconauta belgica Auctorum.] 


1. Indice 100 x - 


N 
NE 50) À =2140 — — 11,96 
(ox 
1 M0 rl 4 
5 N S 
h —7 — = 0,174 RE — 4,7720 
À C) V2 T 
+ : ; 
100 © œ X @i)h (Gi) x A |£ M | © à ne 
— Li i p (Xi): @ (Li Gi — 1 ee due & les TR —— 2 
d o Ver | Ve Via Le 
42 —5 1 0,14 —30,1 —,40 0,02239 0,268 
49 —À4 5 0,71 —23,1 —1,84 0,07362 0,880 
56 —3 17 2,43 —16,1 21698 0,17593 2,104 
63 —? 31 4,43 — 9,1 —0,73 0 30548 3,604 
70 —1 29 4,14 — 2,1 —0,17 0,39281 4 G98 
Ki 0 33 4,71 + 4,9 0,39 0,36958 4,420 
84 1 16 2,29 1169 0,95 0,25403 3,038 
91 À dl 1,57 18,9 1551 0,12766 1,527 
98 Ÿ 0,57 25,9 2,07 0,04098 0,490 
105 4 3 0,43 32,9 2,62 0,01295 0,155 
Écart probable : 0,6745 x 12,54=8,46. 
d à 
2. Indice 100 x -. | 
a 
NEO À — 136,98 se 8,681 
n=A5 RER 
" C3 N 
h 11 0,120 EE — 3,4690 
A C) V2 fe 
2 v 2 
d < CAN À ETC | | 
100 — — %; X; p (@) À © (5) Xi — À d=—— | ——68 9 | —-—20 2 
7 ù V2» oV/2r 
101 +5 2 0,29 —35,98 —2,08 Ô,04598 0,399 
108 —À 3 0,43 —?8,98 —1,68 0,09737 0,845 . 
115 —3 HAS 2157 —?21,98 her 0,17882 1,547 | 
122 —? 17 2,43 —14,98 —0,87 0,27317 rap yal 
129 —1 er 3,86 — 7,98 —0,46 0,35855 3,118 | 
136 0 20 2,86 — 0,98 —0,06 0,39775 3,453 | 
143 1 19 or 6,02 +0,35 0,37483 3,254 
150 è 20 2,86 13,02 0,75 0,30097 2,613 
157 3 12 1,71 20,02 1,16 0,22731 1,973 
164 4 5 0,71 27,02 1,56 0,11836 1,027 
171 5 2 0,29 34,02 1,97 0,05748 0,499 
178 6 4 . 0,57 41,02 2,37 0,02417 0,210 d 
185 7 — — 48,02 2,78 0,00842 0,073 
192 8 — — 55,02 3,18 0,00256 0,022 
199 9 


1 0,14 62:02 3,99 0,00064 0,006 


Écart probable : 0,6754 x 17,28 — 11,55. 
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8. Indice 100 x ©. 
5 (1 

N — 146 À — 170,95 . OP 

n = 13 5 — 22,24 

h — 10 © 26430 PEN a6 

À 5 \E K 
RES x o(æ)h o(%:) 50 Fe Un 1 IN —$ 
9 yat ex alt bre 

108 —(6 3 0,3 —(62,95 000 6,00733 0,048 
118 —5 if 0,1 —52,95 —?,38 0,02358 0,155 
128 —À4 1 0,1 —42,95 —1,93 0,06213 0,408 
138 —3 L 0,1 —32,95 —1,48 0,13356 0,876 
148 —? 23 2,3 —22,95 | 108) 0,23473 1,540 
158 —1 28 2,8 —12,95 = (0:08 0,33699 2,e11 
168 0 2 al — 2,95 —0,13 0,39518 2,592 
178 il 30 3,0 7.05) +0,32 0,37877 2,485 
188 2 19 1,9 +17,05 0,77 0,28292 1,856 
198 3 4 0,7 Sey ADS 1,22 0,18963 1,244 
208 n 5 0,5 +37,05 1,67 0,09911 0,650 
218 6) 4 0,4 +47,05 2,12 0,04228 0,277 
228 6 | 0,3 + 57,05 AN 0,01477 0,097 


Écart probable : 


8 6. ÉTUDE DES ANGLES CARACTÉRISTIQUES. 


Considérant toujours la mesure de la grandeur d comme une bonne estima:. 
tion de l’âge des coquilles, il a été dressé des tables de répartition des couples de 
variables statistiques (d, x) et (d, 8). Ces tables ont révélé une absence de corré- 
lation entre les caractères angulaires 4, G et d aisément explicable et dont la 
signification sera dégagée plus loin. Par contre, les distributions de la fréquence 

des variables 4 et $ sont très révélatrices; elles montrent une concentration nette 
des mesures autour de certaines valeurs (voir Pls, XV à XVIIT et tableaux 


pp. 84 à 86). 
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Gisement de Wérister. 
[Anthraconauta minima Hinp (non LUpWwIG), espèce de référence.] 


1. Angle «. 


N 
N — 421 A — 390,66 — — 67,252 
T 
n = 6 = 60,26 
5 s ë N ee 
DES as 0,150 — — 26,830 
(2 « 
Xi — À 1 Re N ns 
à 0 Xi. Nr) ET ol) le AE e 2 = 6 T0 
Ë F Ver ol 27% 
25 5 3 0,6 —14,66 Dr 0,02596 1,745 
30 = 46 9,2 — 9,66 —1,54 0,12208 8,210 
35 =: 118 23,6 — 4,66 —0,74 0,30323 20,393 
40 0 133 26,6 0,34 +0,05 0,39794 26,762 
45 1 66 13,2 5,34 0,85 0,27789 18,689 
50 2 43 8,6 10,34 1,65 0,10249 6,893 
55 3 10 2,0 15,34 2,45 0,01996 1,342 
60 4 1 0,2 20,34 3,25 0,00205 0,138 
65 5 il 0,2 25,34 4,05 0,00011 0,007 
Écart probable : 0,6745 x60,26=— 40,222. 
2. Angle 6. 
N | 4 
N= 1423 A — 1070,91 — — 45,582 
(ox 
tal o — 90,28 
c N 
À C V2 T 
d2 Le 
LC — 1 { er N TER 
B= 2 X4 o(2% )R p(æi) Gi—A |t—— He 2 = 
g Var ol” 27 
100 —6 2 0,4 Ci —3,00 0,00443 0,202 
105 —5 3 0,6 2201 y; 0,01899 0,866 
110 —4 23 4,6 17,91 É5(108 0,06213 2,832 
115 +) 32 6,4 —12,91 —1,39 0,15189 6,923 
120 =? 57 11,4 201 —0,85 0,27789 12,667 
125 ef 64 12,8 2970 —0,31 0,38008 17,325 
130 0 111 22,2 2,09 0,23 0,38814 17,692 
135 1 67 13,4 7,09 0,76 0,29871 13,616 
140 à 47 9,4 12,09 1,30 0,17137 7,811 
145 3 13 2,6 17,09 1,84 0,07362 3,356 
150 4 4 0,8 22,09 2,38 0,02358 1,075 


Écart probable : 0,6745 x 90,28—60,259. 


_ 
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Gisement d’Angleur. 
[Anthraconauta belgica Auctorum.] 


1. Angle «. 


N 
N — 197 À = 380,37 — — 24,472 
Cy 
h—= 12 os — 80,05 
o N 
h= 5 — — 0,2098 — — 9,76 
A 5 JE 5 
9 Æ it 
2 : 2 
Li — À À SES N LYS 
o(xi) dj — À t < e 2 : ne 2 
k V/2= V2 T 
0,2 —21,37 —2,65 0,01200 0,29 
0,6 —16,37 —?,00 0,05399 1:82 
3,4 —11,37 —1,4 0,14973 3,06 
8,8 — 6,37 —0,7 0,31225 7,64 
12,6 — 1,37 —0,2 0,39104 9,57 | 
jl 7,6 + 3,63 0,4 0,36827 9,01 
2 3,0 + 8,63 Lt 0,21785 5383 
3 1,8 + 13,63 do 0,09405 2,30 
4 1,4 +18,63 2,3 4 0,02833 0,69 
62 D: il 0,2 +23,63 2,9 0,00595 0,15 
67 6 ee > + 28,63 3,6 0,00061 0,01 
72 7 è 0,4 +33,63 4,R 0,00006 = 
Écart probable : 0,6745x80,05=50,43, 
2. Angle $ 
N 
N'— 195 A — 129,062 — — 18,2075 
C} 
= A2 006 
C} N À 
R=5 — = 0,0821 —#,20 
A s V27T 
2 
4 Xi — A 1 pe té N nés 
B— CG o(t)h p(æi) ni ANT — + Re 2 
Ê É V’2x op 27 
103 —5 5 1 —?26,06 —?,5 0,01753 0,32 
108 —5 8 1,6 —21,06 —2,0 0,05399 0,98 
113 —3 5 nl —16,06 —1,5 0,12952 2,35 
118 —2 20 4 —11,06 —1,0 0,24197 4,41 
123 —1 19 3,8 — 6,06 —0,57 0,33877 6,17 
128 0 51 10,2 — 1,06 —0,1 0,39695 7,83 
133 1 38 7,6 3,94 0,4 0,36827 6,71 
138 7) 25 5 8,94 0,8 0,28969 Diet 
143 3 9 1,8 13,94 1,3 0,17137 see 
148 4 9 1,8 18,94 1,8 0,07895 1,44 
153 5 3 0,6 23,94 2,3 0,02833 0,52 
158 6 1 0,2 28,94 LAN 0,01042 0,19 


Écart probable : 0,6745x100,6=70,15 
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Gisement de Stud. 
[Anthraconauta belgica Auctorum.] 


N 
N'— 149 MARS US 5 = 19,426 
n\= 44 s — 60,67 
5 N J 
ES Fe 0,178 = 7,150 
SN 2T 
AR 2 
Xi — A À Ne N san. 
CHE Xi o (æi }h o(æi) %Gi—A |é— 6 "9 |) 
? V 2T os 27 
21 —5 1 0,33 —16,53 —?,48 0,01850 0,359 
24 7 5 1,67 —13,53 —2,03 0,05099 0,990 
on —3 5 1,67 10158 —1,58 0,11464 2,227 
30 —? 18 6 — 7,53 —1,13 0,21075 4,094 
33 —1 8 2,67 — 4,53 —0,68 0,31644 6,147 
26 0 42 14 — 1,53 —0,23 0,38814 7,940 
39 1 35 11,67 1,47 +0,22 0,38911 7,559 
42 2 8 2,67 4,47 0,67 0,31854 6,188 
45 9 15 5 7,47 1,12 . 0,21312 4,140 
48 4 2 0,67 10,47 1,57 0,11650 2,263 
bi 5 8 2:07 13,47 2,08 0,05199 1,010 
54 6 — — 16,47 2,47 0,01899 0,369 
57 fé — en 19,47 2,92 0,00565 0,110 
60 8 2 0,67 22,47 BR 0,00138 0,027 
Écart probable : 0,6745 x 6,67 =4,50. 
2. Angle $. 
N 
NEA A — 1350,99 Pa 14,86 
= A0 & — 90,89 } 
C3 N 
C V2 T 
æi—A| 1 Æ | Len | 
B— œj Xi o(æi)R S (4) %i— A t= 4 6 29 
: : V/or sV/2r 
110 4 ri 0,2 __25,99 263 0,01264 0,188 
115 —4 ti 1,0 —?20,99 —?,12 0,04228 0,268 
120 —3 7 1,4 —15,99 —1,62 0,10755 1,598 
125 —? 15 3,0 —10,99 » —1,11 0,21548 3,208 
130 —1 24 4,8 — 5,99 —0,61 0,33112 4,920 
139 0 28 5,6 —.0;99 —0,10 0,39695 5,899 
140 L 29 5,8 4.01 0,41 0,36665 5,448 
145 e 18 3,6 9,01 0,91 0,26367 3,918 
150 3 16 de 14,01 1,42 0,14569 2,165 
155 A 4 0,8 19,01 1,92 0,06332 0,941 


Écart probable : 


0,6745 x 9,89 = 6,67. 
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F. — CONCLUSIONS BIOLOGIQUES RELATIVES A LA « VARIABILITE 
INDIVIDUELLE » ET AUX « VARIATIONS SPECIFIQUES ». 


Les conclusions biologiques de l'étude biométrique des populations 
d'Anthraconauta des gisements À. — Angleur, À. — Stud et À. — Wérister 
doivent montrer pour chacune d'elles l'existence ou l’absence, d’une part, de 
corrélation entre les caractères mesurés et la « variabilité individuelle » de 
ceux-ci, d'autre part, de « formes variantes », dites spécifiques. A l'instar de 


P. PELSENEER (*), nous considérons la « variabilité morphologique indivi- 


duelle » comme résultant d'un ensemble de caractères qui, étant communs à 
tous les individus, sont donc normaux. Ces formes sont distinctes des « variétés 
spécifiques » portant, elles, des modifications d'ordre plus élevé. Rappelons 
encore que cette variabilité est, comme toujours en paléontologie, celle de la 
coquille, organe extérieur de l'animal, et qu’elle ne reflète pas la nature des 
variations correspondantes de l'organisme. Variabilité individuelle et, éventuel- 
lement, variations spécifiques doivent être examinées en fonction de l’âge (ou 
de la taille) et du point de vue de l'intensité, de l'orientation et de la fréquence. 


La croissance ou « variation-normale » de l’organisme détermine des chan- 
gements de la forme de la coquille, grands-et intenses dans les premières 
phases du développement et qui diminuent de l’état embryonnaire à l’état adulte. 
Il convient donc d'essayer de situer la place qu'occupent les Anthraconauta dans 
la suite ontogénique des formes du groupe auquel elles appartiennent, aux fins 
de retrouver éventuellement leur type morphologique à l’un des stades de cette 
suite. Par là on pourrait retrouver ainsi dans les phases de cette série le sens des 
variations possibles. 

Le développement ontogénique général des Mytilidés, auxquels nous rappor- 
tons les Anthraconauta, a été relativement bien étudié. Le rapprochement du 
genre avec des modiolomorphes ou des mytilimorphes n’a pas, à priori, une 
sionification phyllétique par la similitude morphologique, mais il représente au 
moins un Caractère de convergence très net. Ce phénomène, bien connu, est 
caractéristique, pour de nombreux groupes, du passage de la structure isomyaire 
à la structure anysomyaire chez les formes byssogènes, fixées en position pleuro- 
thétique [R. Avrmony (*)]. L’embryogénie des espèces du groupe des Mytilidés 
a été étudiée par R. J. Jackson (%), F. BernarD (7), R. ANrnony (%), T. C. NeL- 


(54) PELSENEER, P., 1919. 

(55) ANTHONY, R., 1905, pp. 230 et suiv. 
(85) JACKSON, R. T., 1890, pp. 363-364. 
(57) BERNARD, F., 1898, pp. 67-86. 

(58) ANTHONY, R., 1905, pp. 244 et suiv. 
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son (f°) et d’autres encore. Elle a montré, dans le développement de la coquille 
embryonnaire, les phases suivantes, schématisées par la figure 11 : 


Dimm 
C Pen 
; D 
E F 
imm 2mm 


FIG. 11. — Mytilus edulis LINNÉ. 


A. Prodyssoconque définitive. — B. Dyssoconque népionique. 
C. Coquille au stade néologique. — D. Coquille au stade « modioloïide ». 
E et F. Stades mytilloïdes. 


Dès le début de la croissance de la dyssoconque ou coquille vraie, au stade 
népionique, celle-ci développe une asymétrie progressive, passant par une phase 
de modiolisation, pour prendre ensuite le type mytilimorphe. 

Ce processus apportant à la coquille son aspect définitif s'effectue par: 

1° un allongement de la valve du côté postérieur, le côté antérieur s'accrois- 
sant moins vite, sans être cependant frappé d’avortement; 

2° un accroissement du bord cardinal vers l'extrémité postérieure de la 
valve, ou, relativement, déplacement de l’umbo vers l’angle antérieur; 

3° un déplacement angulaire de la position du diamètre de croissance maxi- 
mum (d) par réduction de la valeur angulaire de l'angle «x, ou, relativement, 
allongement oblique de la coquille. 


(5°) NELSON, T. C., 1928, pp. 180-190. 
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Ce processus de différenciation graduelle de la dyssoconque, aux stades 


népionique et néalogique, est propre à différents phyllums d’anysomyaires 


byssogènes. Rencontré dans des groupes divers, il ne signifie pour les uns qu'un 
cas de convergence; pour d’autres, il est le film d’une longue évolution phyllo- 
génétique. Nous sommes enclins à considérer le type anthraconautiforme comme 
le représentant d’un stade de l’ontogénèse des Mytilidés, très voisin de l’état 
définitif auquel ont abouti les Modiola (voir chapitre 1, $ 2). Les variations 
«spécifiques », et à plus forte raison les variations « individuelles », du type 
anthraconautiforme apparaissent être déterminées principalement par la fluc- 
tuation des valeurs angulaires « et $ et de la valeur de l'angle de mytilisation 
[R. ANTHONY ("°)]; corrélativement, elles entraînent des modifications des valeurs 
indicielles. 


Les conclusions de l'étude biométrique des trois populations A. — Angleur, 
A. — Stud, À. — Wérister sont également valables pour chacune d'elles. Les 
voici : 


a) Il ÿ a une bonne corrélation entre le caractère d (diamètre de croissance 
maximum, exprimant la taille, ou l’âge et les caractères a (hauteur), b (lon- 
gueur), g (longueur du bord cardinal). Les valeurs élevées du coefficient de 
corrélation (voir tableau, p. 75) sont voisines de 1. Les valeurs les moins élevées, 
0,864 et 0.886 pour o,, de À. — Angleur et o,, de À. — Stud sont à attribuer à 
une dispersion assez forte provenant de valeurs de a et g éloignées des moyennes. 
On peut donc conclure à l'existence, de relations linéaires entre les caractères 
précités. Les rapports entre éléments de même nom de figures semblables étant 


constants, les indices caractéristiques 100 a 100 F 100 ?.. sont indépendants du 
a 


développement normal de la coquille. De la plus petite taille observée (2 mm.) à la 
plus grande (18,3 mm.), la coquille, au cours du développement, reste semblable 
à elle-même; sa croissance est homothétique dès l’adolescense. H. BoLTon ar 
W. Hp ("), E. Dix et AE. TRuEMman () avaient noté que les petites coquilles 
de À. minima, À. laevis et A. phillipsi, grandes de 2 à 4 mm., paraissent généra- 
lement à contour, plus arrondi et à plus faible obliquité. L'impression de cet état, 
que nous avons eue nous-mêmes, semble, après l'examen biométrique de petites 
formes analogues, n'être qu'une illusion de la vision, caractéristique de l’obser- 
vation de petits segments de droites et de courbes. 

b) L'absence de corrélation entre les caractères angulaires 4, B et le carac- 


tère linéaire d se conçoit aisément, puisque, d’une part (au a), on constate le 
développement homothétique des coquilles au cours de la croissance et que, 


NTHONY, R., 1905, p. 234. 

OLTON, H., 1915,/pp. 676-677. 

IND, W., 1894, p. 272. 

IX, E. et TRUEMAN, A. E., 1931. p. 188. 


CI > 


90 ÉTUDE BIOMÉTRIQUE DES ANTHRACONAUTA 


d'autre part, les distributions de fréquences des variables x et 6 montrent une 


concentration nette des mesures obtenues autour de certaines valeurs. 


c) L'examen du tableau de corrélation et des graphiques de distribution de 
fréquence des indices dénote, dans les limites observées des tailles d'individus, 
comme on pouvait s’y attendre. qu'il n’y a pas de discontinuité de croissance. 
Ceci, compte tenu des stases ou des ralentissements de la croissance dus aux 
phénomènes saisonniers et aux périodes d'activités reproductrices 


d) L'examen du tableau de distribution de fréquence des indices caracté- 
ristiques et des angles «et 6 met en lumière que la « variabilité individuelle », 
toujours existante, présente une dispersion normale, régulière. Les courbes de 
fréquences n’offrant pas d’asymétrie, il n’y a pas de sens observable à cette 
variabilité, pas plus qu'on n’y peut reconnaître une amplitude ou une multi- 
plicité, En un mot, l'allure de la « variabilité individuelle » est normale et il 
n'y a de « variétés spécifiques » dans aucun des trois gisements étudiés. 


e) La dispersion normale des indices prouve que ces individus fossiles 
examinés ont été récoltés dans leur habitat naturel, à l’état grégaire, et n'ont 
pas été rassemblés en gisement, tels certains « mussels band », par action méea- 
nique du milieu. 


G. — COMPARAISON DES INDICES TYPES DES DIFFERENTES 
POPULATIONS POUR LA RECHERCHE D’UNE CONSPECIFICITE 
EVENTUELLE DE LEURS ELEMENTS. 


La comparaison de la spécificité de populations limniques ou dulcicoles, 
appartenant à des gisements différents, appelle le souvenir de cette remarque 
très générale de P. PELSENEER : « La variabilité est inégale entre des individus 
ou groupes d'individus dans une même espèce; elle est la plus grande dans les 
régions où les conditions d'existence sont les plus variées ». Elle est encore 
d'autant plus marquée que l'aire de dispersion de l'espèce est plus grande. C'est 
bien Ià le cas des Anthraconauta du terrain houiller étudiées. Cependant les 
individus des colonies étudiées ici et provenant de trois gisements distincts et 
éloignés présentent une conspécificité très nette. 


L'ajustement des séries statistiques, obtenues suivant la loi de répartition 
de Gauss, s'étant montré très satisfaisant, on a fait usage d’une propriété des 
écarts quadratiques, propre à la loi précitée, et permettant de définir l'intervalle 
quasi certain d'apparition des valeurs. On a admis aussi que la moyenne est la 
valeur la plus probable et que la « variable moyenne » définie plus haut (p. 27) 
obéit à la loi. Cette façon de faire est peut-être plus légitime que d'ajuster la 
série observée suivant la variable statistique elle-même par une répartition de 
Gauss, qui n’est vraie que d’une façon asymptotique. 
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La « variable moyenne » correspondant à un indice caractéristique d’une 
espèce varie de façon quasi certaine dans un intervalle dont on peut déterminer 
l'ampleur. Cet intervalle est d’ailleurs d'autant plus petit que le nombre m de 
mesures est grand; il varie comme _— . Si des mesures analogues sont effectuées 

m 
sur plusieurs indices caractéristiques d’une même espèce, on détermine chaque 


Intervalle quasi-certain de la variable moyenne : A-3V2 &_ — A+3V27 & 
Vm Va. 
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FIG. 12. 


fois l'intervalle certain de leur apparition. Dès lors, si des mesures, portant sur 
d’autres espèces, révèlent des valeurs moyennes situées dans ces intervalles, il 
ést permis de conclure, avec une probabilité d'autant plus grande que m est 
plus grand, à l'identité de ces espèces. 

Les résultats obtenus sur les trois populations étudiées d’A. belgica (Aucto- 
rum) des gisements d’Angleur et de Stud et d'A. minima, du gisement de Wéris- 
ter, espèce de référence, ont été groupés en tableaux relatifs, chaque fois à un 
indice ou à un angle caractéristiques. 
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De façon générale, À et 5 étant la moyenne et la dispersion de l'indice ou 
de l'angle étudié, l'intervalle quasi-certain de la variable moyenne est 


ue (o} Fa (ox 
A 3 V2. RE SAN Den, 
Un Un 
_. Espèce Valeur : : Intervalle 
Indice ouorigine |moyenne| Dispersion | "” | quasi-certain 

A. — Wérister 71,94 10,76 433 de. 69,75 à 74,13 

1005 = A.— Angleur 72,92 14,10 209! de “68,78 à 77,06 
b |A. Stud 72,10 42,54 450 | de 67,75 à 76,45 

d A.— Wérister 135,71 43,38, 433 | de 133,05 à 138,49 

100 X — A.— Angleur 131,26 21,04 209 | de 125,08 à 137,44 
à A.— Stud 136,98 17,28 150 | de 130,99 à 142,97 

b A. — Wérister 179,52 20,51 425 | de 175,30 à 183,74 
100 A.— Angleur 165,27 26,00 201 | de 157,47 à 173,07 
8 A. — Stud 170,95 22,24 146 | de 163,14 à 178,76 

; Espèce Valeur : : . Intervalle 
Angle ou origine | moyenne Dispersion | quasi-certain 

| A.— Wérister 392,66 6° ,26 421 de 382,37 à 40°,95 

œ A.— Angleur 380,37 8°,05 197 | de 352,94 à 400,80 

{ A. — Stud 319,53 6°,67 149 de 35°,20 à 39°,86 
A.— Wérister 4270 ,91 90,28 423 | de 125°,99 à .1290,83 
8 ÀA.— Angleur 1290 ,06 100,60 193 | de 1250,85 à 1320,27 
A,— Stud 1359 ,99 90 ,89 147 | de 1320,52 à 139c°,46 


L'examen comparatif des valeurs moyennes du tableau, transcrites sur le 
graphique de la figure 12, est déjà significatif : la distance maximum entre deux 
de ces valeurs est de l’ordre de grandeur maximum de 8 %; les intervalles quasi- 
certains des « valeurs moyennes » se chevauchent au moins largement. 

IL est remarquable d'arriver à une concordance telle avec des populations 
peu nombreuses, et, bien que l'argument des grands nombres ne puisse être 
invoqué de manière incontestable, on peut conclure avec de très fortes pré- 
somptions à l'identité entre les espèces étudiées, à savoir entre l’Anthraconauta 
belgica Hinn et l’Anthraconauta minima (Hixp, non Lupwic). 


CONCLUSIONS GENERALES. 


Si l'étude biométrique des organismes fossiles ne révèle rien sur la nature 
ou le principe des phénomènes biologiques, « elle permet cependant de résoudre 
des problèmes de spécificité difficilement solubles », et aussi de préciser les 
limites de la « variation spécifique ». 

« Pour que les méthodes statistiques donnent les résultats que nous sommes 
en droit d’en attendre, il est nécessaire qu’un nombre considérable de documents 
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soient accumulés; il est donc désirable que les études sur ce sujet se multiplient. 
Mais, pour que les résultats de ces études soient comparables, il est indispen- 
sable qu'à l'instar des anthropologistes, les paléontologues standardisent leurs 
méthodes » (M. GLIBERT, 1945, p. 4). Au préalable, celles-ci doivent être expo- 
sées, confrontées et expérimentées, ainsi que précisées dans leur signification 
mathématique profonde et dans le sens réel de leurs possibilités pratiques, à la 
fois par le biologiste et le mathématicien, pour justifier et ordonner pleinement 
leur emploi. 

Au surplus, dans de nombreux cas, les grandeurs de l'objet lui-même, qui 
deviennent mesures et sont considérées comme caractères morphologiques, 
demandent à être examinées auparavant dans leur valeur biologique et dans 
leur valeur géométrique. Sur ce dernier point, force nous est de constater que 
les caractères d'une coquille (tels que la longueur, la hauteur, etc.), que les 
zoologistes estiment représentatifs, sont parfois dépourvus du sens précis que 
leur conférerait une définition où une propriété géométrique rigoureuse, Mesurer 
des « longueurs » ou des « largeurs » est une entreprise raisonnable, à condition 
qu on sache définir exactement ces grandeurs. On ne mesure pas de façon bien 
précise une « longueur » de coquille, quand on ne sait pas exactement où elle 
commence ni où elle finit, étant donné que le segment de droite qui la repré- 
sente est limité par des ares. Or, c'est cela qu'on se contente de faire le plus 
souvent. 

A notre avis, à des caractères difficiles à mesurer convenablement, on 
aurait tout intérêt à substituer des caractères définis sans ambiguïté par des 
repères géométriques simples. En nombre suffisant, ceux-ci peuvent permettre, 
par l'application de méthodes statistiques, de déterminer d’abord l'équation de 
la courbe du bord libre de la coquille, puis l'équation de sa surface géométrique. 
Ce dernier problème a déjà été traité par M. Lison (1942, pp. 377 et suiv.). Dans 
un travail ultérieur nous reviendrons sur l'application de ces méthodes qui appa- 
raissent rationnelles. 

Quoi qu'il en soit, même avec les imperfections résultant et de l'état de 
conservation du matériel étudié, et des définitions de certains caractères mesurés, 
l'étude biométrique de l'Anthraconauta belgica Hi et de l'Anthraconauta 
minima Hivo (non Lupwic) a démontré une fois de plus qu'on peut arriver 
à des conclusions à probabilités très élevées qui. pour le biologiste, ont la valeur 
de certitudes. 
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EXPLICATION DE LA PLANCHE A 
A | 


F1G. 1 à 27. — Anthraconauta minima (HINp, non LUupwiG), non Anthraconaula belgic 
Hi. 
Loc. : PI. Chênée, nos 4 et 5, Angleur. 


Fic. 1 à 5 et 6 à 10. — (N°" A 3148, À 230, À 500, À 478, À 454) et (N° À 18, AE 
À 465, A 23, À 22). Suite d'exemplaires représentatifs du type moyen, norma 


Fic. 11 à 16. — (N°° À 348, A 10, A 69, À 458, À 472, À 145). Suite d’exemplan 
présentant une ligne cardinale et une région umbonale bien caractérisées. L 


FiG. 17 à 20. — (N°* À 435, À 453, À 445, À 465). Suite d'exemplaires présentant une 
ornementation bien caractérisée. ‘4 pe 


Fic. 22 à 27. — (N°* A 429, À 451, A 465, À 452, À 439, À 367, À 294). Suite d'exem- 
plaires présentant une zone umbonale du test différenciée. ‘2% 


Phachel Fa un # j 2 


22 (x3) 23 (x3) 24 (x 3) 


{ e® 
E Anthraconauta minima (HiIND, non LUDWIG), non Anthraconauta belgica HIN», 
| Gisements d'Angleur X 
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EXPLICATION DE LA PLANCHE B 


I 


PTS , 
FiG. 1 à 145. — Anthraconauta minima (Hip, non Lupwi@), non Anthraconauta belgica « 
HIND. 


Loc. : Pl. Andenne, n° 19, Stud. 
FiG. 1 à 10. — (N° S 338, S 451, S 35, S 100, S 451, S'493, S 512, S 170, S 29; S20): 
Suite d'exemplaires représentatifs du type moyen, normal. Ù 


Fic. 11 à 15. — (N° S 30, S 233, S 22, S 253, S 222). Suite d'exemplaires présentant 
une ligne cardinale et une région umbonale bien caractérisées. 


Il 


F1G. 16 à 36. — Anthraconauta minima (HIND, non LUDWIG). 


Loc.: Charbonnage de Wérister; siège Josée, étage de 242 m., bouveau Nord, toit de 
première veinette sous Grosse, niveau n° 43, 


FiG. 16 à 21 et 22 à 27. — (N°° J 82; J 258, J 242, J 189, J 244) et (N° J 406 06/0 
J 367, J 225, J 11). Suite d'exemplaires représentatifs du type moyen, normal. M 


F1G. 18 et 29. — (N° J 486, J 184). Suite d'exemplaires présentant une ornementation 
bien caractérisée. 


Fi1G. 30, 31 et 36. — (N°° J 82, J 407, J 200). Suite d'exemplaires présentant une zone 
umbonale du test différenciée. 


FiG. 32 à 35. — (N° J 89, J 477, J 267, J 419). Suite d'exemplaires présentant une , 
ligne cardinale ef une région umbonale bien caractérisées. 


Pléiche, 5./| ”AHAMNES 


I. - Anthraconauta minima (HIND non LUDWIG), non Anfhraconauta belgica HIND. Gisement de Stud, - 


Il, - Anthraconauta minima (HiND non LUDWIG). Gisement de Wérister. 
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